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Abstract 

Using AdS/CFT, we study the addition of an arbitrary number of backreacting flavors to 
the Klebanov-Witten theory, making many checks of consistency between our new Type IIB 
plus branes solution and expectations from field theory. We study generalizations of our 
method for adding flavors to all M = 1 SCFTs that can be realized on D3-branes at the 
tip of a Calabi-Yau cone. Also, general guidelines suitable for the addition of massive flavor 
branes are developed. 
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1 Introduction 



The AdS/CFT conjecture originally proposed by Maldacena pQ and refined in [21 [5J is one 
of the most powerful analytic tools for studying strong coupling effects in gauge theories. 
There are many examples that go beyond the initially conjectured duality and first steps 
in generalizing it to non-conformal models were taken in [4]. Later, very interesting de- 
velopments led to the construction of the gauge-string duality in phenomenologically more 
relevant theories i.e. minimally or non-supersymmetric gauge theories [5j. 

Conceptually, a clear setup for duals to theories with few SUSY's is obtained by breaking 
conformality and (perhaps partially) supersymmetry, deforming M = 4 SYM with relevant 
operators or VEV's. The models put forward in [5] are very good examples of this. 

Even when there are important technical differences, in the same line of thought, we can 
consider the model(s) developed by Klebanov and a distinguished list of physicists: Witten 
[6], Nekrasov [7], Tseytlin [8], Strassler [9], Herzog and Gubser [10] and Dymarsky and 
Seiberg [TTJ. In these papers (and many extensions of them), a far reaching idea has been 
developed, namely to flow to a confining field theory with minimal SUSY starting from an 
J\f = 1 SCFT with a product gauge group SU(N C ) x SU(N C ), bifundamental chiral matter 
and a quartic superpotential for the chiral superfields0 The superconformal field theory 
described above rules the low energy dynamics of N c D3-branes at the tip of the conifold. 
Then conformality is broken by the addition of fractional branes, that effectively unbalance 
the ranks of the gauge groups [TJ |8] . A "duality cascade" starts and the flow to the IR leaves 
us with a confining field theory [9]. Subtleties related to the last steps of the cascade have 
been discussed in pUl 12] • All this interesting physics is very nicely described with great 
detail in [15] . 

In this paper we will concentrate on a nonconformal theory without cascade. The starting 
point is a Type IIB solution dual to an SU(N C ) x SU(N C ) Af = 1 SCFT also known as the 
Klebanov- Witten field theory /geometry. One of the aims of the paper is to add an arbitrary 
large number of flavors to each of the gauge groups. The addition of fundamental degrees of 
freedom is an important step toward the understanding of QCD-like dynamics, in different 
regions of the space of parameters. 

A very fructiferous idea used to add flavors to different field theories (using the string dual) 
was described in [TJ] and then applied to various backgrounds, 'flavoring' different dual field 
theories, in many subsequent publications (for a complete list see citations to [14J). As it 
was clearly stated in the original paper, the procedure spelled out in [TJ] consists in the 
addition of a finite number Nf of spacetime filling flavor D7-branes to the N c — > oo color 
D3-branes extending in the Minkowski directions, and the usual decoupling limit (g s — > 0, 
N c — *■ oo, g s N c fixed) of the D3-branes is performed, keeping the number Nf of flavor branes 

1 It is obvious that such a field theory is non-renormalizable and must be thought of as the IR of some 
UV well defined theory. In [5] a UV completion in terms of an orbifolded Af = 2 field theory is given. 
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fixed. Then the D3-branes generate the geometry and the flavor branes only minimize 
their worldvolume Dirac-Born-Infeld action in this background without deforming it. This 
is the probe limit. In the dual description they are considering the addition of a finite 
number of flavors Nf to the large N c gauge theory, in the strict double scaling 't Hooft limit 
(gYM —> 0, N c — > oo, A = gy M N c fixed). In the lattice literature this is called the 'quenched' 
approximation: the dynamics of the colors and its effect on the flavors is completely taken 
into account, but the backreaction of the flavors onto the colors is neglected. In the probe 
limit this approximation becomes exact. 

It is interesting to go beyond this 'quenched' or 'non-backreacting' approximation and see 
what happens when one adds a large number of flavors, of the same order of the number 
of colors, and the backreaction effects of the flavor branes are considered. Indeed, many 
phenomena that cannot be captured by the quenched approximation, might be apparent 
when a string backreacted background is found. 

In this paper we will propose a Type IIB dual to the field theory of Klebanov and Witten, 
in the case in which a large number of flavors (Nf ~ N c ) is added to each gauge group. 
We will also present interesting generalizations of this to cases describing different duals to 
M = 1 SCFT's constructed from D3-branes placed at singularities. 

Let us briefly describe the procedure we will follow, inspired mostly by the papers [15, 
[T6l [T7] and more recently [HI [19]. In those papers (dealing with the addition of many 
fundamentals in the non-critical string and Type IIB string respectively), flavors are added 
into the dynamics of the dual background via the introduction of Nf spacetime filling flavor 
branes, whose dynamics is given by a Dirac-Born-Infeld action. This dynamics is intertwined 
with the usual Einstein-like action of IIB and a new solution is found, up to technical 
subtleties described below. 

1.1 Generalities of the Procedure Used 

To illustrate the way flavor branes will be added, let us start by considering the background 
of Type IIB supergravity that is conjectured to be dual to the Klebanov- Witten field theory: 
an M = 1 SCFT with gauge group SU (N c ) x SU (N c ), two chiral multiplets of bifundamental 
matter Ai, B i} i — 1,2 and a (classically irrelevant) quartic superpotential 

W = XTi(A i B j A k Bi) e ik e jl . (1.1) 
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The dual Type IIB background reads 

(2 2 

ds 2 = h{r)- 1/2 dx\ 3 + h(r) 1/2 \ dr 2 + T —Y + sin 2 0i dip 2 ) + ^-(# + V cos^ dtp- 

6 z — ' 9 ^— ' 

k i=l,2 i=l,2 

F 5 = -(l + ^AAdlifr)- 1 
277T(7 s iV ca /2 

= 4^ (L2) 

with constant dilaton and all the other fields in Type IIB supergravity vanishing. The set of 
coordinates that will be used in the rest of the paper is x M = {x°, x 1 , x 2 , x 3 , r, 0\, <f\, 9 2 , <p 2 , ip}- 
For the sake of brevity, in the following we will take units is which g s = 1, a' = 1. 

We will add Nj spacetime filling D7-branes to this geometry, in a way that preserves some 
amount of supersymmetry. This problem was studied in pOj |21] for the conformal case and 
in [22] for the cascading theory. These authors found calibrated embeddings of D7-branes 
which preserve (at least some fraction of) the supersymmetry of the background. We will 
choose to put two sets of D7-branes on the surfaces parametrized by 

= {x°, x , x 2 , x 3 , r, 6*2, y>2, $i = const. tpi = const. , 

£2 = {%°, % , x 2 , x 3 , r, 0i, (fx, ip} 9 2 = const. ip 2 = const. . (1.3) 

Note that these two configurations are mutually supersymmetric with the background. More- 
over, since the two embeddings are noncompact, the gauge theory supported on the D7's has 
vanishing 4d effective coupling on the Minkowski directions; therefore the gauge symmetry 
on them is seen as a flavor symmetry by the 4d gauge theory of interest. The two sets of 
flavor branes introduce a U(Nf) x U(Nf) symmetry^, the expected flavor symmetry with 
massless flavors. The configuration with two sets (two branches) can be deformed to a single 
set, shifted from the origin, that represents massive flavors, and realizes the explicit breaking 
of the flavor symmetry to the diagonal vector-like U(Nf). Our configuration (eq. (11.31) ) for 
probes is nothing else than the z\ — holomorphic embedding of [2T] . 

We will then write an action for a system consisting of type IIB supergravitjjjjl plus D7- 
branes described by their Dirac-Born-Infeld action (in Einstein frame): 



5 = -V / A^fe-^M^-^i^i 2 -!^! 2 

ZK-lq J VI Z 4 

N f N f 



+ 



;i-4) 



2 The diagonal axial U(l) is anomalous 

3 The problems with writing an action for Type IIB that includes the self-duality condition are well known. 
Here, we just mean a Lagrangian from which the equations of motion of Type IIB supergravity are derived. 
The self-duality condition is imposed on the solutions. 
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where we have chosen the normalization |-F p | 2 = p \F p F p (G~ 1 ) p . 

Notice that we did not excite the worldvolume gauge fields, but this is a freedom of the 
approach we adopted. Otherwise one may need to find new suitable ^-symmetric embeddings. 

These two sets of D7-branes are localized in their two transverse directions, hence the 
equations of motion derived from (11.41) will be quite complicated to solve, due to the presence 
of source terms (Dirac delta functions). 

But we can take some advantage of the fact that we are adding lots of flavors. Indeed, 
since we will have many (Nf ~ iV c ) flavor branes, we might think about distributing them 
in a homogeneous way on their respective transverse directions. This 'smearing procedure' 
boils down to approximating 



N f 



T 7 N 



N r 



All 

T 7 N f 
Air 



d w x e 4, sin 0,- 



-G 



(0 



Volpi) + Vol(Y 2 ) 



AC 8 



with Vol{Yj) = sin6*j dQ{ A dipi the volume form of the S 2, s. 
This effectively generates a ten dimensional action 



2k| J L 2 



+ 



Mi 

47T 



L — <■ v 4^ 

i=l,2 



Vol(Y x ) + Vol(Y 2 



AC* 



;i.5) 



;i.6) 



We can derive in the smeared case the following (not so involved) equations of motion, 
coming from the action (11.61) : 



Rmn — -GmnR 



D M d 



AM 



dF b 



i (d M (fid N (j) - ^G MN d P <; 



1 „2</,/p(l) p (l) l r |p(l)|2V 
2 e y M N ~ 2 G ^I^ I ) + 



e 2<f >\F l \ 2 + 




(5) p (5)PQR5 
PQRS N 

2k 2 T 7 N 



96 Mp Q RS N 



MN 



e* 1 > sin^ 

J^G 4vr ^ 

v i=l,2 



-2k 2 10 T 7 
. 



4tt 



Volpi) + Vol{Y 2 ) 



-G 



W 



;i-7) 



The modified Bianchi identity is obtained through Fi = — e 2 ^ * Fg, and comes from the WZ 
part of the action ( 11.61) . The contribution to the stress-energy tensor coming from the two 
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sets of Nf D7 flavor branes is given by 

/,M - "~ - -Tf-^>.-^-\/-^"^"" w . (i.> 



:.\/.\ _ 2 'no 5S flavor _ N f < \ - u>(i ; , .,/_/„.i'i/.'.'>-/ A -u A -.v 

\J —G SGmn 

where a, /3 are coordinate indices on the D7. In the subsequent sections we will solve the 
equations of motion (ll.7jl - fll.8D and will propose that this Type IIB background is dual to the 
Klebanov-Witten field theory when two sets of Nf flavors are added for each gauge group. 
We will actually find BPS equations for the purely bosonic background, by imposing that 
the variations of the dilatino and gravitino vanish. We will verify that these BPS first-order 
equations solve all the equations of motion (11.71) . 

Let us add some remarks on some important points about the resolution of the system. 
First of all, it is clear from the Bianchi identity of F\ in (11.7ft that we will not be able to 
define the axion field Cq on open subsets. 

Regarding the solution of the equations of motion, we will proceed by proposing a deformed 
background ansatz of the form 

ds 2 = h~ 1/2 dxl 3 + h l/2 \dr 2 + — ^{d9 2 + sin 2 0*6,$) + — (# + Yl C0S M<ft) 2 ) 

^ ^ i=l,2 ^ i=l,2 ' 

F 5 = (l + *)d 4 x AKdr (1.9) 
N f 

F\ = - — (dip + cos 9 2 dif2 + cos 9\ dipi) . 

Thanks to the smearing procedure, all the unknown function h, f, g, K and the dilaton 
only depend on the radial coordinate r. 

The Bianchi identity for the five-form field-strength gives 

Kh 2 e A9+s = 27nN c , (1.10) 

and we will obtain solutions to (11.71) by imposing that the BPS equations derived from the 
vanishing of the gravitino and gaugino variations and the Bianchi identities are satisfied. 
These will produce ordinary first-order equations for f(r), g(r), h(r), K(r), <p(r). We will 
also be able to derive these BPS equations from a superpotential in the reduction of Type 
IIB supergravity. 

We will study in detail the dual field theory to the supergravity solutions mentioned above, 
making a considerable number of matchings. The field theories turn out to have positive 
/5-function along the flow, exhibiting a Landau pole in the UV. In the IR we still have a 
strongly coupled field theory, which is "almost conformal" . We will also generalize all these 
results to the interesting case of a large class of different Af = 1 SCFTs, deformed by the 
addition of flavors. In particular we will be able to add flavors to every gauge theory whose 
dual is AdS$ x M 5 , where M 5 is a five-dimensional Sasaki-Einstein manifold. New solutions 
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will be found that describe the 'unflavored' case, making contact with old results. Finally, a 
possible way of handling the massive flavor case is undertaken. 

We have explained the strategy we adopt to add flavors, so this is perhaps a good place 
to discuss some interesting issues. The reader might be wondering about the 'smearing 
procedure' discussed above, what is its significance and effect on the dual gauge theory, 
among other questions. It is clear that we smear the flavor branes just to be able to write a 
10-dimensional action that will produce ordinary (in contrast to partial) differential equations 
without Dirac delta functions source terms. 

The results we will show and the experience obtained in [T7JEIE] show that many properties 
of the flavored field theory are still well captured by the solutions obtained following the 
procedure described above. It is not clear what important phenomena on the gauge theory 
we are losing in smearing, but see below for an important subtlety. 

One relevant point to discuss is related to global symmetries. Let us go back to the weak 
coupling (g s N c — > 0) limit, in which we have branes living on a spacetime that is the product 
of four Minkowski directions and the conifold. When all the flavor branes of the two separate 
stacks ( j 1 .31) are on top of each other, the gauge symmetry on the D7's worldvolume is given 
by the product U(Nf) x U (Nf). When we take the decoupling limit for the D3-branes a' — > 0, 
with fixed g s N c and keeping constant the energies of the excitations on the branes, we are 
left with a solution of Type IIB supergravity that we propose is dual to the Klebanov-Witten 
field theory with Nf flavors for both gauge groups [21]. In this case the flavor symmetry is 
U(Nf) x U(Nf), where the axial U(l) is anomalous. This background would be for sure very 
involved, since it would depend on the coordinates (r, 61,62), if the embeddings of the two 
stacks of D7-branes are 6\ = and 62 = 0, respectively. When we smear the Nf D7-branes, 
we are breaking U(Nf) — > U(l) Nf (see Figured]). 

There is one important point to contrast with [T7j. In that paper, a smearing is also 
proposed but it is argued that the dual field theory (a version of Af = 1 SQCD with a 
quartic superpotential in the quark superfields) possesses U(Nf) global flavor symmetry. As 
in all backgrounds constructed on wrapped branes, the effects of the Kaluza-Klein modes 
play an important role and the dual field theory behaves as 4-dimensional only in the far 
IRH In this regime, when the internal manifold shrinks to small size, for energies below 
this inverse size we do not effectively see the breaking U(Nf) — > U(l) Nf . In contrast, the 
backgrounds obtained by placing D-branes at conical singularities, like [6]- [10] as well as our 
solution, describe a four dimensional field theory all along the flow. 

It might be interesting for the reader to note that the papers in the line of [H] are working 
in the context of 't Hooft expansion [21]. When the ratio Nf/N c is very small, one can ignore 
the backreaction effects of the flavor branes on the geometry. This is the dual version to 
the suppression of effects that include the running of fundamentals in internal loops. Even 

4 For a detailed study of the role and dynamics of the KK modes in wrapped brane setups, see [23] , 
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Figure 1: We see on the left side the two stacks of Nj flavor-branes localized on each of their 
respective S 2, s (they wrap the other S 2 ). The flavor group is clearly U(Nf) x U(Nf). After 
the smearing on the right side of the figure, this global symmetry is broken to U(l) N f 1 x 
Uil)^- 1 x U(1) B x U(1) A . 

when these fundamentals are massless, their effects while running in loops are suppressed 
by a factor of 0(Nf/N c ). But in the strict 't Hooft limit, if the number of flavors is kept 
fixed, the corrections due to the quantum dynamics of quarks exactly vanish. In the cases 
considered in [T71 ITS] , the ratio above is of order one and we are working on the so called 
Veneziano's topological expansion [25J. New physics (beyond the 't Hooft limit) is captured 
by Veneziano's proposal: we will be able to see this in the present paper that considers the 
backreaction of the flavor branes, just like was observed in [T7J [TB|, in contrast with the 
papers that worked in the 't Hooft approximation as proposed in |14j . 

Another point that is worth elaborating on is whether there is a limit on the number of 
D7-branes that can be added. Indeed, since a D7-brane is a co dimension-two object (like a 
vortex in 2 + 1 dimensions) its gravity solution will generate a deficit angle; having many 
seven branes, will basically "eat-up" the transverse space. This led to the conclusion that 
solutions that can be globally extended cannot have more than a maximum number of twelve 
D7-branes [26] (and exactly twenty-four in compact spaces). In this paper we are adding 
a number Nf —>■ oo of D7-branes, certainly larger that the bound mentioned above. Like 
in the papers [271 12H], we will adopt the attitude of analyzing the behavior of our solutions 
and we will see that they give sensible results. But there is more than that: the smearing 
procedure distributes the D7's all over this 2-dimensional compact space, in such a way that 
the equation for the axion-dilaton is not the one in the vacuum at any point. This avoids 
the constraint on the number of D7-branes, which came from solving the equation of motion 
for the axion-dilaton outside sources. 

Finally, we must emphasize that this is not the first paper that deals with the D3/D7 
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system in the context of "AdS/CFT with flavors". Indeed, very good papers have been 
written where this problem was faced looking for a solution where the flavor branes are 
replaced by fluxes in terms of the Type IIB supergravity fields, dilaton and an axion (0, Co). 
The BPS equations for the D3/D7 system in cases preserving 8 supercharges were written 
in [271 EE], a partially explicit solution of the equations of motion in the presence of sources 
was found in [29] for the orbifold case, more interesting geometrical aspects were discussed 
in [30] and an involved solution was found in [31J, where some matching with gauge theory 
behavior was attempted^ 

The papers [27j EH EH E0J, El] were written with the idea of letting the flavor branes 
backreact. One qualitative difference with respect to what we explained above is that the 
authors of [271 Ell [301 EE] consider the case in which D3-branes are added in the background 
produced by D7-branes and solve the Laplace equation, in this case for the deformation 
introduced by the D3's. In contrast, we consider the background produced by the D3-branes 
and we deform it to take into account the "smeared" backreaction of the D7-branes. The 
two procedures are different. 

One advantage of the approach proposed in [17J is that the flavor degrees of freedom 
explicitly appear in the DBI action that allows the introduction of SU(Nf) gauge fields in 
the bulk that are dual to the global symmetry in the dual field theory, while it is difficult to 
see how they will appear in a Type IIB solution that only includes RR fluxes. Our approach 
produces a simple SUSY solution to ( 11. 7ft and the analysis of gauge theory effects is simple to 
do. Besides, the proposal of [17] used in the present work is the natural continuation of the 
many successful results obtained in papers in the line of [TS] . Indeed, we are just following 
the idea of [H] for a large number of flavor branes. 



1.2 Organization of This Paper 

This paper is organized in two main parts. In part I we will present the addition of flavors to 
the Klebanov-Witten solution. A detailed analysis of the supergravity plus branes solutions 
and the study of the dual gauge theory, as read from the above mentioned solutions, is 
performed. A reader mainly interested with the line of research, but who does not want to 
go in full details, should be happy reading this introduction, part I and Appendix O 

The readers who intend to work on this subject and want to study these results in more 
technical detail or want to appreciate the beauty and generality in our formalism are referred 
to part II. Also in part II the reader will find a sketch of how to deal with massive flavors 
using these techniques. 

Those readers who are not attracted by the physics of flavor using AdS/CFT techniques, 
but just want to learn about some new solutions (born out of our 'deformed backgrounds' 

5 Even though slightly unrelated to the D3/D7 system, we cannot resist here to mention the beautiful 
solution found by Cherkis and Hashimoto for a localized D2/D6 system [32] . 
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as described above), should read the introduction and the appendix B. 
Some other appendices complement nicely our presentation. 

The section of conclusions includes also a summary of results and proposes future directions 
that the interested reader might want to pursue. 

2 Part I: Adding Flavors to the Klebanov-Witten Field 
Theory 

2.1 What to Expect from Field Theory Considerations 

In this first part we will address in detail the problem of adding a large number of backreact- 
ing non-compact D7-branes to the Klebanov-Witten Type IIB supergravity solution, which 
describes D3-branes at the tip of the conifold. Before presenting the solution and describing 
how it is obtained, we would like to have a look at the dual field theory, and sketch which 
are the features we expect. 

For this purpose, we consider the case of probe D7-branes, and mainly summarize what 
was pointed out in |21j. The conifold is a non-compact Calabi-Yau 3-fold, defined by one 
equation in C 4 : 

-21-22 - -^3-24 = . (2.1) 

Since this equation is invariant under a real rescaling of the variables, the conifold is a 
real cone, whose base is the Sasaki-Einstein space T 1 ' 1 |6l [33]. It can be shown that T 1 ' 1 
is a U(l) bundle over the Kahler-Einstein space S 2 x S 2 , and that its isometry group is 
SU{2) x SU(2) x U(l). 

Klebanov and Witten [6] obtained an interesting example of gauge/gravity duality by 
placing a stack of N c D3-branes at the apex of the conifold. The branes source the RR 5- 
form flux and warp the geometry, giving the Type IIB supergravity solution (11. 2p . The dual 
field theory, describing the IR dynamics on the worldvolume of the branes, has gauge group 
SU(N C ) x SU(N C ) and matter fields A^, Bi, i = 1,2 which transform in the bifundamental 
representations (N c , N c ) and (N c , N c ) respectively. The theory has also a quartic superpo- 
tential Wkw = A Tr(AiBjAkBi) e lk e^ 1 . The field theory is Af = 1 superconformal, and the 
anomaly-free U(l) R-symmetry of the superconformal algebra is dual to the U(l) isometry 
of the fiber in T 1,1 , generated by the so-called Reeb vector. In the algebraic definition (12.11) 
it is realized as a common phase rotation of the four coordinates: Zi — > e~ ta Zi. 

The addition of flavors, transforming in the fundamental and antifundamental representa- 
tions of the gauge groups, can be addressed by including probe D7-branes into the geometry, 
following the procedure proposed in [TJ]. This was done in |21j . where the embedding of the 
flavor branes and the corresponding superpotential for the fundamental and antifundamen- 
tal superfields were found. The D7-branes have four Minkowski directions parallel to the 
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stack of D3-branes transverse to the conifold, whereas the other four directions are embed- 
ded holomorphically in the conifold. In particular, D7-branes describing massless flavors can 
be introduced by considering the holomorphic noncompact embedding z\ — 0. The flavors, 
which correspond to 3-7 and 7-3 strings, are massless because the D7-branes intersect the 
D3-branes. Note that the D7-branes have two branches, described by z\ = z% = and 
z\ = = 0, each one corresponding to a stack. The presence of two branches is required 
by RR tadpole cancellation: in the field theory this amounts to adding flavors in vector-like 
representations to each gauge group, hence preventing gauge anomalies. The fundamental 
and antifundamental chiral superfields of the two gauge groups will be denoted as q, q and 
Q, Q respectively, and the gauge invariant and flavor invariant superpotential proposed in 
[2TJ is 

W = W KW + W f , (2.2) 

where 

W KW = A TriAiBtAjBi) e^e kl (2.3) 

is the SU(2) x SU(2) invariant Klebanov-Witten superpotential for the bifundamental fields. 
For a stack of flavor branes, it is conventional to take the coupling between bifundamentals 
and quarks at a given point of S 2 as 

W f = h ~q a A x Q a + h 2 Q a B iga . (2.4) 

This coupling between bifundamental fields and the fundamental and antifundamental flavors 
arises from the D7 embedding z\ = 0. The explicit indices are flavor indices. This superpo- 
tential, as well as the holomorphic embedding z\ = 0, explicitly breaks the SU(2) x SU(2) 
global symmetry (this global symmetry will be recovered after the smearing). 

The field content and the relevant gauge and flavor symmetries of the theory are summa- 
rized in Table [1] and depicted in the quiver diagram in Figure [21 
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3/4 


1 


1 


Q 


(N c , 1) 


(hN f ) 


(1,1) 


3/4 


-1 


-1 


Q 


(l,iV c ) 


(hN f ) 


(1,1) 


3/4 


1 





Q 


(1,N C ) 


(N f A) 


(1,1) 


3/4 


-1 






Table 1: Field content and symmetries of the KW field theory with massless flavors. 

The U(1)r R-symmetry is preserved at the classical level by the inclusion of D7-branes 
embedded in such a way to describe massless flavors, as can be seen from the fact that the 
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Figure 2: Quiver diagram of the Klebanov-Witten gauge theory with flavors. Circles are 
gauge groups while squares are non-dynamical flavor groups. 



equation Z\ = is invariant under the rotation Z{ — > e~ ta Zi and the D7 wrap the R-symmetry 
circle. Nevertheless the U(1)r turns out to be anomalous after the addition of flavors, due to 
the nontrivial Cq gauge potential sourced by the D7. The baryonic symmetry U(1)b inside 
the flavor group is anomaly free, being vector-like. 

As was noted in [21], the theory including D7-brane probes is also invariant under a 
rescaling Zi — > (3zi, therefore the field theory is scale invariant in the probe approximation. 
In this limit the scaling dimension of the bifundamental fields is 3/4 and the one of the flavor 
fields is 9/8, as required by power counting in the superpotential. Then the beta function 
for the holomorphic gauge couplings in the Wilsonian scheme is 



= = > ft = <^A? . (2.5) 

with Aj = gfN c the 't Hooft couplings. In the strict planar 't Hooft limit (zero order in 
Nf/N c ), the field theory has a fixed point specified by the afore- mentioned choice of scaling 
dimensions, because the beta functions of the superpotential couplings and the 't Hooft 
couplings are zero. As soon as Nf/N c corrections are taken into account, the field theory has 
no fixed points for nontrivial values of all couplings. Rather it displays a "near conformal 
point" with vanishing beta functions for the superpotential couplings, but non-vanishing 
beta functions for the 't Hooft couplings. In a Nf/N c expansion, formula (12.51) holds at 
order Nf/N c if the anomalous dimensions of the bifundamental fields A,- and Bj do not get 
corrections at this order. A priori it is difficult to expect such a behavior from string theory, 
since the energy-momentum tensor of the flavor branes will induce backreaction effects on 
the geometry at linear order in Nf/N c , differently from the fluxes, which will backreact at 
order (N f /N c ) 2 . 
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Moreover, since we are adding flavors to a conformal theory, we can naively expect a 
Landau pole to appear in the UV. Conversely, we expect the theory to be slightly away from 
conformality in the far IR. 

2.2 The Setup and the BPS Equations 

The starting point for adding backreacting branes to a given background is the identification 
of the supersymmetric embeddings in that background, that is the analysis of probe branes. 
In [20], by imposing /t-symmetry on the brane world- volume, the following supersymmetric 
embeddings for D7-branes on the Klebanov-Witten background were found: 

£° = {x°, x, x 2 , x 3 , r, 9 2 , ip%, ip} 6>i = const. <p\ = const. (2 6) 

^2 = {x°,x 1 ,x 2 ,x 3 ,r,9i,Lfi,-ip} 9 2 = const. <£> 2 = const. 

They are precisely the two branches of the supersymmetric embedding z\ — first proposed 
in [21]. Each branch realizes a U(Nf) symmetry group, giving the total flavor symmetry 
group U(Nf) x U(Nf) of massless flavors (a diagonal axial U(1)a is anomalous in field theory, 
which is dual to the corresponding gauge field getting massive in string theory through Green- 
Schwarz mechanism). We choose these embeddings because of the following properties: 
they reach the tip of the cone and intersect the color D3-branes; wrap the U(1)r circle 
corresponding to rotations ip — > ip + a; are invariant under radial rescalings. So they realize 
in field theory massless flavors, without breaking explicitly the U(1)r and the conformal 
symmetry. Actually, they are both broken by quantum effects. Moreover the configuration 
does not break the Z 2 symmetry of the conifold solution which corresponds to exchanging 
the two gauge groups. 

The fact that we must include both the branches is due to D7-charge tadpole cancellation, 
which is dual to the absence of gauge anomalies in field theory. An example of a (non- 
singular) 2-submanifold in the conifold geometry is T>2 = {9\ = # 2 , ipi = 2tt — ip2, ip = 
const, r = const}. The charge distributions of the two branches are 

UJ w = y 5^(e 1 ,<pi)de 1 Adip 1 uo {2) = Y ^(#2,^2)^2 a , (2.7) 

— 'TV/ — <N f 

where the sum is over the various D7-branes, possibly localized at different point s, and a cor- 
rectly normalized scalar delta function (localized on an 8-submanifold) is 5^ (x) \f— ~G%/ y/—G. 
Integrating the two D7-charges on the 2-submanifold we get: 

f J 1 ) = -N f [ lu {2) = N f . (2.8) 

Thus, whilst the two branches have separately non- vanishing tadpole, putting an equal num- 
ber of them on the two sides the total D7-charge cancels. This remains valid for all (non- 
singular) 2-submanifolds. 
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The embedding can be deformed into a single D7 that only reaches a minimum radius, 
and realizes a merging of the two branches. This corresponds to giving mass to flavors and 
explicitly breaking the flavor symmetry to SU(Nf) and the R-symmetry completely. These 
embeddings were also found in |20j. 

Each embedding preserves the same four supercharges, irrespectively to where the branes 
are located on the two 2-spheres parameterized by (6*i, cpx) and (9 2 , (p 2 )- Thus we can smear 
the distribution and still preserve the same amount of supersymmetry. The 2-form charge 
distribution is readily obtained to be the same as the volume forms on the two 2-spheres in 
the geometry, and through the modified Bianchi identity it sources the flux We expect 
to obtain a solution where all the functions have only radial dependence. Moreover we were 
careful in never breaking the Z 2 symmetry that exchanges the two spheres. The natural 
ansatz is: 

ds 2 = h(r)- l/2 dxi fi + h(r) 1/2 |dr 2 + 

e 2g(r) g2/(r) ^ (2.9) 

+ — — 2^ (d9 2 + sin 2 6i dip 2 ) + — — (dip + ^ cos^ d^) 2 \ 

i=l,2 '" i=l,2 ^ 

<p = <p(r) (2.10) 
F 5 = K{r) h(r) 3 '* (e^^ - e 6 ^ 102 ^ (2.11) 

F 1 = - l (d^ + cos 0id^i + cos 6 2 dip 2 ) = ^ h(r)- 1/4 e - f{r) (2.12) 

47T 4:71 

dF l = L (sin 6 X dQ x A dip x + sin 6 2 dQ 2 A <kp 2 ) , (2.13) 

47T ' 

where the unknown functions are h(r), g(r), f(r), <ft(r) and K(r). The angular coordinates 
8i are defined in [0, tt] while the others have fundamental domain ipi e [0, 2tt) and ip G [0, 4tt) 

6 The modified Bianchi identity of Fx is obtained from the Wess-Zumino action term with Fx = — e^ 2 ^*Fg. 
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with appropriate patching rules^. The vielbein is: 

e* = hT 1 ' A dx i e r = h l ' 4 dr 

= 7t hl/4e9 d9i ^ = ^ /4 e 9 sin0 f <fy>i (2.14) 

e i> _ h llA e f {dip + cos 61 d<pi + cos 9 2 d<p 2 ) . 

With this ansatz the field equation d(e 2< ^ * Fi) = is automatically satisfied, as well as the 
self-duality condition F 5 = *F 5 . The Bianchi identity dF 5 = gives: 

Kh 2 e A9+f = 27nN c , (2.15) 

and K{r) can be solved. The previous normalization comes from Dirac quantization of the 
D3-brane charge: 

/ F 5 = 24,T 3 N c = {2n) 4 N c , (2.16) 
Jt 1 - 1 



using a suitable orientation for the volume form of the T 1 ' 1 space and the fact that Vo^T 1 ^ 



16^3 
27 



We impose that the ansatz preserves the same four supersymmetries as the probe D7- 
branes on the Klebanov-Witten solution. With this purpose, let us write the supersymmetric 
variations of the dilatino and gravitino in type IIB supergravity. For a background of the 
type we are analyzing, these variations are: 



6 e \ = ~r M (d M <l>-ie+F£ [ y 
S e iP M = V M e + z jF$ e + ^F^ RST T^ RST T M e , (2.17) 



where we have adopted the formalism in which e is a complex Weyl spinor of fixed ten- 
dimensional chirality (see Appendix IA"]) . It turns out (see Section 13.21) that the Killing 
spinors e (which solve the equations 5 e A = 5 e ipM = 0) in the frame basis (I2.14p can be 
written as: 

e = h~* 77 (2.18) 



7 The correct patching rules on T 1 ' 1 in the coordinates of p.2p are: 

In fact the space is a t/(l) fibration over S 2 x S 2 . The first identification is just the one of the fiber. On 
the base 2-spheres we must identify the angular variables according to ipi = tp^ + 2-7T, but this could be 
accompanied by a shift in the fiber. To understand it, draw the very short (in proper length) path around 
the point 6% — 0: 9% <C 1, ipi = t = 4tt — if> with t £ [0, 27r] a parameter along the path. To make it closed, a 
rotation in <p\ must be accompanied by an half-rotation in ip. This gives the second identification. 
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where r\ is a constant spinor which satisfies 



T^i^s 77 = -if] 



rV^i 77 = r^^a r] = ir] , T r ^7] = -it] . (2.19) 

Moreover, from (12. 1 7[) we get the following system of first-order BPS differential equations: 



/' = e -/(3-2e 2/ - 29 ) - 



3iV 



57T 



3iV 



47T 



(2.20) 



k /i' = -27nN c e' f - Aa 



Notice that taking Nf = in the BPS system (I2.20p we simply get equations for a 
deformation of the Klebanov-Witten solution without any addition of flavor branes. Solving 
the system we find both the original KW background and the solution for D3-branes at a 
conifold singularity, as well as other solutions which correspond on the gauge theory side to 
giving VEV to dimension 6 operators. These solutions were considered in [341 135] . and are 
shown to follow from our system in appendix [Bl 

In order to be sure that the BPS equations (I2.20j) capture the correct dynamics, we have 
to check that the Einstein, Maxwell and dilaton equations are solved. This can be done even 
before finding actual solutions of the BPS system. We checked that the first-order system 
(I2.20p (and the Bianchi identity) in fact implies the second order Einstein, Maxwell and 
dilaton differential equations. An analytic general proof will be given in Section 13.31 In the 
coordinate basis the stress-energy tensor ( 11. 81) is computed to be: 



3iV 



2tt 



T, 



N 



f 



<Pi<Pi 



24vr 



-2g 



4e 2/ cos 2 6i + 3e 2g sin 2 & 



T 

1 r . 



3iV 



T, 



2tt 



To 



N 



■in 



T, 



Iplp 



.^l e 4>+2f-2g cosdi cos02 

6n 

.^^+2/-2 Scos0 . 
Nf_ e<t>+ 2f-2g _ 

Qir 



(2.21) 



It is correctly linear in Nf. We did not explicitly check the Dirac-Born-Infeld equations for 
the D7-brane distribution. We expect them to be solved because of K-symmetry (supersym- 
metry) on their world-volume. 



Solution with General Couplings 

We can generalize our set of solutions by switching on non- vanishing VEVs for the bulk gauge 
potentials C2 and B 2 . We show that this can be done without modifying the previous set 
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of equations, and the two parameters are present for every solution of them. The condition 
is that the gauge potentials are flat, that is with vanishing field-strength. They correspond 
thus to (higher rank) Wilson lines for the corresponding bundles. 
Let us switch on the following fields: 

C 2 = cu 2 B 2 = bu 2 , (2.22) 

where the 2-form uj 2 is Poincare dual to the 2-cycle T> 2 : 

F> 2 = {#i = 2 , if\ = 2n — if 2 , i> = const, r = const} (2.23) 

uj 2 = — ( sin 9 1 ddi A dipi — sin 9 2 dd 2 A d<p 2 ) , [ uj 2 = 1 . (2.24) 

8vr v ' J V2 

We see that = and = 0. So the supersymmetry variations are not modified, neither 
are the gauge invariant field-strength definitions. In particular the BPS system (12.201) does 
not change. 

Consider the effects on the action (the argument is valid both for localized and smeared 
branes). It can be written as a bulk term plus the D7-brane terms: 

S = S Mk -T 7 J d 8 £e^-det(G 8 + F)+T 7 J J ^ Cq A e*] g , (2.25) 

with T — B 2 + 2ita' F is the D7 gauge invariant field-strength, and hat means pulled-back 
quantities. To get solutions of the K-symmetry conditions and of the equations of motion, 
we must take F such that 

T = B 2 + 2Tra'F = . (2.26) 

Notice that there is a solution for F because B 2 is flat: dB 2 = dB 2 = 0. With this choice 
K-symmetry is preserved as before, since it depends on the combination T. The dilaton 
equation is fulfilled. The Bianchi identities and the bulk field-strength equations of motion 
are not modified, since the WZ term only sources Cg. The energy momentum tensor is not 
modified, so the Einstein equations are fulfilled. The last steps are the equations of B 2 and 
A\ (the gauge potential on the D7). For this notice that they can be written: 

^ = 2to '^ = (2.27) 



5B 2 5B 2 

The first is solved by T = since in the equation all the terms are linear or higher order 
in T . This is because the brane action does not contain terms linear in J 7 , and this is true 
provided Cq = (which in turn is possible only if C 2 is flat). The second equation then 
reduces to Ss 5 b ^ k = 0, which amounts to d(e~^ * H 3 ) = and is solved. 
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As we will see in Section 12. 5[ being able to switch on arbitrary constant values c and b 
for the (flat) gauge potentials, we can freely tune the two gauge couplings (actually the two 
renormalization invariant scales A's) and the two theta angles [6], [36]. This turns out to 
break the Z2 symmetry that exchanges the two gauge groups, even if the breaking is mild 
and only affects C2 and B2, while the metric and all the field-strength continue to have that 
symmetry. However this does not modify the behavior of the gauge theory. 

2.3 The Solution in Type IIB Supergravity 

The BPS system (12.201) can be solved through the change of radial variable 

e~ f dr = dp. (2.29) 



j d _ d 
dr dp 



We get the new system: 



g = e 2f ~ 29 (2.30) 

/ = 3 - 2e 2 ^ - (2.31) 

4,= 3 -§Le* (2.32) 

h = -27tt N c e~ 49 , (2.33) 

where derivatives are taken with respect to p. 

Equation (12.321) can be solved first. By absorbing an integration constant in a shift of the 
radial coordinate p, we get 

! ' => p<0. (2.34) 



3N f p 

The solution is thus defined only up to a maximal radius Pmax = where the dilaton 
diverges. As we will see, it corresponds to a Landau pole in the ultraviolet (UV) of the 
gauge theory. On the contrary for p — > — 00, which corresponds in the gauge theory to the 
infrared (IR), the string coupling goes to zero. Note however that the solution could stop at 
a finite negative pmin due to integration constants or, for example, more dynamically, due 
to the presence of massive flavors. Then define 

u = 2f-2g => w = 6(l-e u ) + -, (2.35) 

P 

whose solution is 

e« - ~ 6pe6P (2 36) 
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The constant of integration c\ cannot be reabsorbed, and according to its value the solution 
dramatically changes in the IR. A systematic analysis of the various behaviors is presented 
in Section 12.41 The value of c\ determines whether there is a (negative) minimum value for 
the radial coordinate p. The requirement that the function e u be positive defines three cases: 

-1 < ci < -> pmin < P < 
ci = — > — oo < p < 
ci > — > -oo < p < . 

In the case — 1 < c\ < 0, the minimum value pmin is given by an implicit equation. It can 
be useful to plot this value as a function of C\\ 

Pmin 



= (1 - 6pmin) e 6 " MIN + d 



As it is clear from the graph, as C\ — > — 1 + the range of the solution in p between the IR 
and the UV Landau pole shrinks to zero size, while in the limit C\ — * 0~ we no longer have 
a minimum radius. 

The functions g(p) and f(p) can be analytically integrated, while the warp factor h(p) 
and the original radial coordinate r(p) cannot (in the particular case c\ = we found an 
explicit expression for the warp factor). By absorbing an irrelevant integration constant into 
a rescaling of r and a;°' 1 ' 2 ' 3 ; we get: 

1/6 

(2.37) 



1 - 6p)e 6p + a 
6p~e 3p 



1 - 6p)e 6p + ci 



-1/3 



Hp) 

r(p) 



27itN, I , 
o 

p 

J 



-Mi 



+ c 2 



(2.38) 
(2.39) 

(2.40) 



This solution is a very important result of our paper. We accomplished in finding a super- 
gravity solution describing a (large) Nf number of backreacting D7-branes, smeared on the 
background produced by D3-branes at the tip of a conifold geometry. 

The constant c\ and c 2 correspond in field theory to switching on VEV's for relevant 
operators, as we will see in Section 12.5.31 Moreover, in the new radial coordinate p, the 
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metric reads 

i i f e 2fl ~ 2/ x - 1 1 

^ ^ i=l,2 ^ i=l,2 J 

(2.41) 

2.4 Analysis of the Solution: Asymptotics and Singularities 

We perform here a systematic analysis of the possible solutions of the BPS system, and 
study the asymptotics in the IR and in the UV. In this section we will make use of the 
following formula for the Ricci scalar curvature, which can be obtained for solutions of the 
BPS system: 

■a v, . ._ _ r -a V, 1 

(2.42) 



An 



2.4.1 The Solution with c x = 



7 _|_ j 3N f c 2fl-2/+(fr 

An 



Although the warp factor h(p) cannot be analytically integrated in general, it can be if the 
integration constant C\ is equal to 0. Indeed, introducing the incomplete gamma function, 
defined as follows: 



we can integrate 



%) = -277riV c y dp z +c 2 = 



= lnN c (^T[ 1 -,- 2 - + A P ] + c 2 ^ (2-44) 
27 

~ — nN c (-6p)- 2/3 e- 4p for p -> -oo . 

The warp factor diverges for p — > — oo, and the integration constant c 2 disappears in the IR. 
Moreover, if we integrate the proper line element ds from a finite point to p = — oo, we see 
that the throat has an infinite invariant length. 

The function r(p) cannot be given as an analytic integral, but using the asymptotic be- 
havior of e-f for p — > — oo we can approximately integrate it: 

r(p) ~ 6 1 / 6 [(-p) 1 ^ + lr[i -p]] + c 3 (2.45) 
in the IR. Fixing r — > when p — > — oo we set c 3 = 0. We approximate further on 

r(p) ~ (-6p) 1/6 e p . (2.46) 
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Substituting r in the asymptotic behavior of the functions appearing in the metric, we find 
that up to logarithmic corrections of relative order 1/| log(r)|: 



e 9(r) ^ g/(r) ^ r 



h{r) 



27ttN c 1 (2-47) 



4 r 4 

Therefore the geometry approaches AdS^ x T 1,1 with logarithmic corrections in the IR limit 
p — > — oo. 



2.4.2 UV Limit 



The solutions with backreacting flavors have a Landau pole in the ultraviolet (p — > 0~), 
since the dilaton diverges (see (12.341) ). The asymptotic behaviors of the functions appearing 
in the metric are: 



,29 



,2/ 



h 



;i+ci 

-6p(l 



,1/3 



6p 2 



1 + Ci 



o(p 3 



Cl 



-2/3 



f 6p + C(p 
c 2 + 277riV c (l + c 1 )- 2/3 [-p- 



1+Ci 



P 3 + <V 



(2.48) 
(2.49) 
(2.50) 



Note that we have used (I2.39P for the warp factor. One concludes that h(p) is monotonically 
decreasing with p; if it is positive at some radius, then it is positive down to the IR. If the 
integration constant c 2 is larger than zero, h is always positive and approaches c 2 at the 
Landau pole (UV). If c 2 = 0, then h goes to zero at the pole. If c 2 is negative, then the warp 
factor vanishes at Pmax < before reaching the pole (and the curvature diverges there). 
The physically relevant solutions seem to have c 2 > 0. 

The curvature invariants, evaluated in string frame, diverge when p — > 0~, indicating 
that the supergravity description cannot be trusted in the UV. For instance the Ricci scalar 
R ~ (— p)~ 5 / 2 if c 2 7^ 0, whereas R ~ (— p)~ 3 if c 2 = 0. If c 2 < 0, then the Ricci scalar 
R ~ (Pmax - py 1/2 when p -> p MAX . 



2.4.3 IR Limit 

The IR (p — > — oo) limit of the geometry of the flavored solutions is independent of the 
number of flavors, if we neglect logarithmic corrections to the leading term. Indeed, at the 
leading order, flavors decouple from the theory in the IR (see the discussion below eq. (12. 5p ). 
The counterpart in our supergravity plus branes solution is evident when we look at the BPS 
system (12.201) : when p — ► — oo the term disappears from the system, together with all 
the backreaction effects of the D7-branes (see Appendix [C] for a detailed analysis of this 
phenomena), therefore the system reduces to the unflavored one. 
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ci = 

The asymptotics of the functions appearing in the metric in the IR limit p —>■ — oo are: 

e 9 ~ e 7 ~ (-6p) 1/6 e p (2.51) 
97 

ft ~ tL^ c ^ p y^ e -^ . (2.52) 

Formula (I2.42p implies that the scalar curvature in string frame vanishes in the IR 
limit: R {s) ~ (-p)~ 1/2 -> 0. An analo gous but lengthier formula for the square of the 
Ricci tensor gives 

R mnR (S)MN = ^^- c (-P) + <W - oo, (2.53) 

thus the supergravity description presents a singularity and some care is needed when 
computing observables from it. The same quantities in Einstein frame have limiting 
behavior ~ (-p)~ 1/2 -> and R { ^ N R^ MN -> 640/(277ri\Q. 

ci > 

The asymptotics in the limit p — > — oo are: 

e 9 ~ q /6 (2.54) 
~ c x 1/3 (-Qp) 1/2 e 3p (2.55) 
/i ~ 27nN c Ci 2/3 (-p) . (2.56) 

Although the radial coordinate ranges down to — oo, the throat has a finite invariant 
length. The Ricci scalar in string frame is R ~ (— p)~ 3 e~ 6p — > — oo. 

ci < 

In this case the IR limit is p — > pmin- The asymptotics in this limit are: 

e" ~ ( - 6pMiNe 6/,MIN ) V6 (6p - 6 Pm in) 1/6 (2.57) 

' ~ ( - 6p M iNe 6pMIN ) 1/6 (6p - 6 Pm in)" 1/3 (2.58) 
h ~ const. > . (2.59) 



e 



The throat has a finite invariant length. The Ricci scalar in string frame is R ~ 
(p - Pmin)~ 1/3 -> oo. 

Using the criterion in [37], that proposes the IR singularity to be physically acceptable if g tt 
is bounded near the IR problematic point, we observe that these singular geometries are all 
acceptable. Gauge theory physics can be read from these supergravity backgrounds. We call 
them "good singularities". 
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2.5 Detailed Study of the Dual Field Theory 



In this section we are going to undertake a detailed analysis of the dual gauge theory features, 
reproduced by the supergravity solution. The first issue we want to address is what is the 
effect of the smearing on the gauge theory dual. 

As we wrote above, the addition to the supergravity solution of one stack of localized non- 
compact D7-branes at Z\ — put in the field theory flavors coupled through a superpotential 
term 

W = X TriABkAjBt) e ij e kl + h x q a A x Q a + h 2 Q a B iga , (2.60) 

where we explicitly wrote the flavor indices a. For this particular embedding the two branches 
are localized, say, at Q\ = and 9 2 = respectively on the two spheres. One can exhibit a 
lot of features in common with the supergravity plus D7-branes solution: 

• the theory has U(Nf) x U(Nf) flavor symmetry (the diagonal axial U(1)a is anomalous), 
each group corresponding to one branch of D7's; 

• putting only one branch there are gauge anomalies in QFT and a tadpole in SUGRA, 
while for two branches they cancel; 

• adding a mass term for the fundamentals the flavor symmetry is broken to the diagonal 
U(Nf), while in SUGRA there are embeddings moved away from the origin for which 
the two branches merge. 

The SU (2) x SU (2) part of the isometry group of the background without D7's is broken 
by the presence of localized branes. It amounts to separate rotations of the two S 2 in 
the geometry and shifts the location of the branches. Its action is realized through the 
superpotential, and exploiting its action we can obtain the superpotential for D7-branes 
localized in other places. The two bifundamental doublets Aj and Bj transform as spinors 
of the respective SU(2). So the flavor superpotential term for a configuration in which the 
two branches are located at x and y on the two spheres can be obtained by identifying two 
rotations that bring the north pole to x and y. There is of course at/(l)x U(l) ambiguity in 
this. Then we have to act with the corresponding 577(2) matrices U x and U y on the vectors 
(Ai, A 2 ) and (Bi, B 2 ) (which transform in the (2, 1) and (1, 2) representations) respectively, 
and select the first vector component. In summary we can write § 

W f = h q x [U x ] Q x + h 2 Q y [U y ^ ] q y , (2.61) 

where the notation q x , Q x stands for the flavors coming from a first D7 branch being at x, 
and the same for a second D7 branch at y. 



8 In case the two gauge couplings and theta angles are equal, we could appeal to the Z2 symmetry that 
exchanges them to argue \h\\ = \h%\, but no more because of the ambiguities. 
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To understand the fate of the two phase ambiguities in the couplings hi and h 2} we appeal 
to symmetries. The U(l) action which gives (q, q, Q, Q) charges (1, — 1, —1, 1) is a symmetry 
explicitly broken by the flavor superpotential. The freedom of redefining the flavor fields 
acting with this U(l) can be exploited to reduce to the case in which the phase of the two 
holomorphic couplings is the same. The U(l) action with charges (1, 1, 1, 1) is anomalous 
with equal anomalies for both the gauge groups, and it can be used to absorb the phase 
ambiguity into a shift of the sum of Yang-Mills theta angles Q\ M + Q\ M (while the difference 
holds steady). This is what happens for D7-branes on flat spacetime. The ambiguity we 
mentioned amounts to rotations of the transverse K 2 space, whose only effect is a shift of Cq. 
As we show in the next section, the value of Cq is our way of measuring the sum of theta 
angles through probe D(-l)-branes. Notice that if we put in our setup many separate stacks 
of D7's, all their superpotential £7(1) ambiguities can be reabsorbed in a single shift of Cq. 

From a physical point of view, the smearing corresponds to put the D7-branes at different 
points on the two spheres, distributing each branch on one of the 2-spheres. This is done 
homogeneously so that there is one D7 at every point of S 2 . The non-anomalous flavor 
symmetry is broken from U(1)b x SU(Nf)n x SU(Nf)i (localized configuration) to U(1)b X 
U(l)y f 1 x U(l)^ f 1 (smeared configuration)^] 

Let us introduce a pair of flavor indices (x, y) that naturally live on S 2 x S 2 and specify the 
D7. The superpotential for the whole system of smeared D7-branes is just the sum (actually 
an integral) over the indices (x, y) of the previous contributions: 



q y 



W = \ Tr(AiB k AjBi) e ij e kl + h x j ^ d 2 x q x [u x f^M ] Q x + h 2 J d 2 y Q y [u y g 

(2.62) 

Again, all the U(l) ambiguities have been reabsorbed in field redefinitions and a global shift 

of e\ M + e\ M . 

In this expression the SU(2)a x SU(2)b symmetry is manifest: rotations of the bulk fields 
Aj, Bj leave the superpotential invariant because they can be reabsorbed in rotations of the 
dummy indices (x, y). In fact, the action of SU{2)a x SU(2) b on the flavors is a subgroup 
of the broken U(Nf) x U(Nf) flavor symmetry. In the smeared configuration, there is a D7- 
brane at each point of the spheres and the group SU (2) 2 rotates all the D7's in a rigid way, 
moving each D7 where another was. So it is a flavor transformation contained in U(Nf) 2 . By 
combining this action with a rotation of Ai and Bi, we get precisely the claimed symmetry. 

Even if written in an involved fashion, the superpotential f 1 2 . 6 2 j) does not spoil the features 
of the gauge theory. In particular, the addition of a flavor mass term still would give rise to 



9 The axial U (1) which gives charges (1, 1,-1, —1) to one set of fields (q X7 q x , Q x , Q x ) coming from a single 
D7, is an anomalous symmetry. For every D7-brane we consider, the anomaly amounts to a shift of the same 
two theta angles of the gauge theory. So we can combine this U(l) with an axial rotation of all the flavor 
fields, and get an anomaly free symmetry. In total, from Nf D7's we can find Nf — 1 such anomaly free axial 
U(l) symmetries. 
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the symmetry breaking pattern 



U(1) B x U{l)y f - X x £7(1) 



k JV>-l 



£7(1) 



A'/' 



2.5.1 Holomorphic Gauge Couplings and /5-functions 

In order to extract information on the gauge theory from the supergravity solution, we need 
to know the holographic relations between the gauge couplings, the theta angles and the 
supergravity fields. These formulae can be properly derived only in the orbifold IR 1,3 x C x 
C 2 /Z2, where string theory can be quantized, by considering fractional branes placed at 
the singularity. The near-horizon geometry describing the IR dynamics on a stack of N 
regular branes at the singularity is AdS 5 x 5 ,5 /Z 2 . The dual gauge theory is an Af = 2 
SU(N) x SU(N) SCFT with bifundamental hypermultiplets. In Af = 1 language, an TV = 2 
vector multiplet decomposes into a vector multiplet and a chiral multiplet in the adjoint of the 
gauge group, whereas a bifundamental hypermultiplet decomposes into two bifundamental 
chiral multiplets. Klebanov and Witten [6] recognized that giving equal (but opposite) 
complex mass parameters to the adjoint chiral superfields of this Af = 2 SCFT, an RG 
flow starts whose IR fixed point is described by the gauge theory dual to the AdS 5 x T 1 ' 1 
geometry. 

In the Af = 2 orbifold theory, the holographic relations can be derived exactly. The result 
is the following: 
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B 2 — tt (mod 27r) 
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C 2 (mod27r) 
C 2 (mod 2tt) 



(2.63) 
(2.64) 
(2.65) 
(2.66) 



where the integrals are performed over the 2-sphere that shrinks at the orbifold fixed point 
and could be blown-upo The ambiguity in (12.641) is the 2n periodicity of f s B 2 which 
comes from the quantization condition on H3 (if fractional branes are absent). A shift of 
2n amounts to move to a dual description of the gauge theory0 The ambiguities of RR 
fields are more subtle: the periodicities in (12.651) and (12.661) correspond to the two kinds of 



10 Actually, we would find the opposite sign in the Co term in the formulas (|2.65p and (|2.66|) for the 6 
angles. We are not sure about that sign. At any rate, with this minus sign the R-anomaly computation of 
the supergravity backgrounds of [5T] and of this paper match exactly the field theory computations. 

In the KW theory, this is Seiberg duality. Notice that the periodicity must fail once flavor fields are 
added. 
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Figure 3: Unit cell of the lattice of Yang-Mills 9 angles and RR fields integrals. 



fractional D(-l)-branes appearing in the theory. The angles 9{ and Q\ come from the 
imaginary parts of the action of the two kinds of fractional Euclidean D(-l) branes. Both of 
them are then defined modulo 2tt in the quantum field theory: 



iYM 



YM\ 



qYM 



+ 2tt) . 



(2.67) 



On the string theory side the periodicities exactly match: an Euclidean fractional D(-l)- 
brane enters the functional integral with a term exp{ — ^ + Hence the imaginary 

part in the exponent is defined modulo 2ir in the quantum string theory. The identification 
(I2.67P of the field theory translates on the string side in: 



2 ) = (irCo + 7T , — 
S 2 Z7C J S 2 



, J- / C 3 +*) = (icC + ir, J- f Q 
Ztt J S 2 Z-k J S 2 



(2.68) 



The lattice is shown in figure [H The vectors of the unit cell drawn in the figure are the ones 
defined by fractional branes. 

From figure [3] and ( 12.681) we can see that: 



ttCq = vrCn + 27T 



(2.69) 



This is indeed the identification that arises from considering a regular D(-l) brane, which 
can be seen as a linear superposition of the two kinds of fractional D(-l)-branes. Notice that 
the closed string field Cq in this orbifold has periodicity 2, differently from the periodicity 1 
in flat space. This is due to the fact that in the orbifold the fundamental physical objects 
are the fractional branes. 



12 We have written the complexified gauge coupling instead of the supergravity fields for the sake of brevity: 
the use of the dictionary is understood. 
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Usually in the literature the afore-mentioned holographic relations were assumed to hold 
also in the conifold case. Strassler remarked in [13] that for the conifold theory the formulae 
for the sum of the gauge couplings and the sum of theta angles need to be corrected. We 
expect that the formula for the sum of theta angles is correct as far as anomalies are con- 
cerned, since anomalies do not change in RG flows. Instead the formula (I2.63H may need 
to be corrected in the KW theory: in general the dilaton could be identified with some 
combination of the gauge and superpotential couplings. 

Let us now make contact with our supergravity solution. In the smeared solution, since 
dF\ 7^ at every point, it is not possible to define a scalar potential Co such that Fi = dC$. 
We by-pass this problem by restricting our attention to the non-compact 4-cycle defined by 
{p, if), 6*i = 9 2 ,<pi = 2n — ip 2 } [2] (note that it wraps the R-symmetry direction ip), so that 
we can pull-back on it and write 

F[ ss = ^ # (2.70) 

4:71 



and therefore 



Now we can identify: 



Co e// = ^(^-^o). (2.71) 



f^e-^-^p (2.72) 

^H yJI = -fW-W, (2-73) 

where we suppose for simplicity the two gauge couplings to be equal (g± — g 2 = g). The 
generalization to an arbitrary constant B 2 is straightforward since the difference of the inverse 
squared gauge couplings does not run. Although, as discussed above, one cannot be sure of 
the validity of (I2.72j) . we can try to extract some information. 

Let us first compute the /3-function of the gauge couplings. The identification (I2.63j) allows 
us to define a "radial" /3-function that we can directly compute from supergravity [38J: 

& s 1*4 = 4± = . (2.74) 

dp g 2 dp 4 

(Compare this result with eq. (j2.5p ). The physical /3-function defined in the field theory is 
of course: 

% = jJ^K > ( 2 - 75 ) 
ir ologf g 2 

where p is the subtraction scale and A is a renormalization group invariant scale. In order 
to get the precise field theory /^-function from the supergravity computation one needs the 
energy-radius relation p = p(^), from which (3 = (3^ dp/dlog^. In general, for non- 
conformal duals, the radius-energy relation depends on the phenomenon one is interested in 
and accounts for the scheme-dependence in the field theory. 
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Even without knowing the radius-energy relation, there is some physical information that 
we can extract from the radial /5-function ( 12.741) . In particular, being the energy-radius 
relation p = p(j) monotonically increasing, the signs of the two beta functions coincide. 

In our case, using r = 4 and eq. f)2.46p . one gets matching between (I2.5P and (I2.74p . 

2.5.2 R-symmetry Anomaly and Vacua 

Now we move to the computation of the U(1)r anomaly. On the field theory side we follow 
the convention that the R-charge of the superspace Grassmann coordinates is R\d\ = 1. This 
fixes the R-charge of the gauginos R[X] = 1. Let us consider an infinitesimal R-symmetry 
transformation and calculate the U{1)r — SU (N c ) — SU (N c ) triangle anomaly. The anomaly 
coefficient in front of the instanton density of a gauge group is RfT[R,^]^ where the sum 
runs over the fermions /, Rf is the R-charge of the fermion and T[TZ^] is the Dynkin index 
of the gauge group representation TZ^ the fermion belongs to, normalized as T[TZ^ wnd ''] = 1 
and T[7^ a *' ) ] = 2N C . Consequently the anomaly relation in our theory is the following: 

M = 3^ (J^*T + G^GT) , (2.76) 

or in other words, under a U(1)r transformation of parameter e, for both gauge groups the 
theta angles transform as 

OjM 6 ym _ N± e (2 J7) 

On the string/gravity side a U(1)r transformation of parameter e is realized (in our 
conventions) by the shift ip — > ip + 2e. This can be derived from the transformation of the 
complex variables (12.11) . which under a U(1)r rotation get zi — ► e t£ Zi, or directly by the 
decomposition of the lOd spinor e into 4d and 6d factors and the identification of the 4d 
supercharge with the 4d spinor. By means of the dictionary ( 12. 73ft we obtain: 

qYM + Q YM _^ B YM + B YM (2.78) 

in perfect agreement with ( 12. 77ft . 

The U(1)r anomaly is responsible for the breaking of the symmetry group, but usually 
a discrete subgroup survives. Disjoint physically equivalent vacua, not connected by other 
continuous symmetries, can be distinguished thanks to the formation of domain walls among 
them, whose tension could also be measured. We want to read the discrete symmetry 
subgroup of U(1)r and the number of vacua both from field theory and supergravity. In 
field theory the U(1)r action has an extended periodicity (range of inequivalent parameters) 
£ G [0, 8ir) instead of the usual 2ti periodicity, because the minimal charge is 1/4. Let us 
remark however that when s is a multiple of 2ir the transformation is not an R-symmetry, 
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since it commutes with supersymmetry. The global symmetry group contains the baryonic 
symmetry U(1)b as well, whose parameter we call a G [0, 27r), and the two actions U(1)r 
and U(1)b satisfy the following relation: Wj?(47r) = Wb(vt). Therefore the group manifold 
U(1)r x U(1)b is parameterized by e G [0, 47r), a G [0, 2iz) (this parameterization realizes 
a nontrivial torus) and U(1)b is a true symmetry of the theory. The theta angle shift 
( I2.77P allows us to conclude that the U(1)r anomaly breaks the symmetry according to 
U(1) r xU(1)b — > Z Nf xU(l) B , where the latter is given by e = Arnr/Nf (n = 0, 1, . . . ,Nf—l), 
a G [0,2tt). 

Coming to the string side, the solution for the metric, the dilaton and the field strengths 
is invariant under arbitrary shifts of if). But the nontrivial profile of Co, which can be probed 
by D(-l)-branes for instance, breaks this symmetry. The presence of DBI actions in the 
functional integral tells us that the RR potentials are quantized, in particular Co is defined 
modulo integers. Taking the formula ( 12. 71 ft and using the periodicity An of ip, we conclude 
that the true invariance of the solution is indeed Z Nf . 

One can be interested in computing the domain wall tension in the field theory by means 
of its dual description in terms of a D5-brane with 3 directions wrapped on a 3-sphere (see 
[39] for a review in the conifold geometry). It is easy to see that, as in Klebanov-Witten 
theory, this object is stable only at r = (p —>■ — oo), where the domain wall is tensionless. 

2.5.3 The UV and IR Behaviors 

The supergravity solution allows us to extract the IR dynamics of the KW field theory with 
massless flavors. Really what we obtained is a class of solutions, parameterized by two 
integration constants c\ and Ci- Momentarily, we will say something about their meaning 
but anyway some properties are independent of them. 

The fact that the /3-function is always positive, with the only critical point at vanishing 
gauge coupling, tell us that the theory is irreparably driven to that point, unless the super- 
gravity approximation breaks down before {c.\ < 0), for instance because of the presence of 
curvature singularities. Using the p coordinate this is clear-cut. In cases where the string 
coupling falls to zero in the IR, the gravitational coupling of the D7 to the bulk fields also 
goes to zero and the branes tend to decouple. The signature of this is in equation (12.311) of 
the BPS system: the quantity e^Nf can be thought of as the effective size of the flavor back- 
reaction which indeed vanishes in the far IR. The upshot is that flavors can be considered 
as an "irrelevant deformation" of the AdS$ x T 1 ' 1 geometry. 

The usual technique for studying deformations of an AdS$ geometry is through the GKPW 
[21, [3] formula in AdS/CFT. Looking at the asymptotic behavior of fields in the AdS 5 effective 
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theory!^! 



5$ = ar A ~ 4 + cr- A , 



(2.79) 



we read, on the CFT side, that the deformation is H = He ft + clO with c = (O) the VEV 
of the operator corresponding to the field $, and A the quantum dimension of the operator 
O. Alternatively, one can compute the effective 5d action and look for the masses of the 
fields, from which the dimension is extracted with the formula: 

A = 2 + V4 + m 2 , (2.80) 

with the mass expressed in units of inverse AdS radius. We computed the 5d effective action 
for the particular deformations e^ r \ e 9 ^ and <p(r) and including the D7-brane action terms 
(the details are in Section [3]). After diagonalization of the effective Kahler potential, we got 
a scalar potential V containing a lot of information. First of all, minima of V correspond to 
the AdS^ geometries, that is conformal points in field theory. The only minimum is formally 
at = 0, and has the AdS$ x T 1,1 geometry. Then, expanding the potential at quadratic 
order the masses of the fields can be read; from here we deduce that we have operators of 
dimension 6 and 8 taking VEV, and a marginally irrelevant operator inserted^ 

The operators taking VEV where already identified in [HI [35] . The dimension 8 operator is 
TrF 4 and represents the deformation from the conformal KW solution to the non-conformal 
3-brane solution. The dimension 6 operator is a combination of the operators Tr(W Q VV Q ) 2 
and represents a relative metric deformation between the S 2 x S 2 base and the U(l) fiber 
of T 1 ' 1 . The marginally irrelevant insertion is the flavor superpotential, which would be 
marginal at the hypothetic AdS§ (conformal) point with = 0, but is in fact irrelevant 
driving the gauge coupling to zero in the IR and to very large values in the UV. Let us add 
that the scalar potential V can be derived from a superpotential W, from which in turn the 
BPS system (I2.20p can be obtained. 

Since in the IR the flavor branes undergo a sort of decoupling, the relevant deformations 
dominate and their treatment is much the same as for the unflavored Klebanov-Witten 
solution [HI [35] [13]. We are not going to repeat it here, and we will concentrate on the case 
Ci = C2 = 0. The supergravity solution flows in the IR to the AdS§ x T 1 ' 1 solution (with 
corrections of relative order 1/| log(r) |). On one hand the R-charges and the anomalous 

13 Notice that usually the GKPW prescription or the holographic renormalization methods are used when 
we may have flows starting from a conformal point in the UV. In this case, our conformal point is in the IR 
and one may doubt about the validity in this unconventional case. See Section 6 in the paper [40] for an 
indication that applying the prescription in an IR point makes sense, even when the UV geometry is very 
far away from AdS^ x M5. We thank Kostas Skenderis for correspondence on this issue. 

14 To distinguish between a VEV and an insertion we have to appeal to the first criterium described in eq. 
(jiTrgi) and below. 
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dimensions tend to the almost conformal values: 
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(2.81) 



1q,Q 
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Using the formula for the /^-function of a superpotential dimensionless coupling: 



/% = H-3+E # ( 1 +t) 



(2.82) 



where $ are the fields appearing in the superpotential term, we obtain that the total su- 
perpotential (12.621) is indeed marginal. On the other hand the gauge coupling flows to zero. 
Being at an almost conformal point, we can derive the radius-energy relation through rescal- 
ings of the radial and Minkowski direction, getting r = p/A. Then the supergravity beta 
function coincides with the exact (perturbative) holomorphic /3-function (in the Wilsonian 



If we are allowed to trust the orbifold relation (I2.63P relating gauge coupling constants 
and dilaton, we conclude that the gauge coupling flows to zero in the IR. This fact could 
perhaps explain the divergence of the curvature invariants in string frame [4], as revealed by 
(I2.53p . The field theory would enter the perturbative regime at this point. However, it is 
hard to understand why the anomalous dimensions of the fields are large while the theory 
seems to become perturbative. For this reason, we question the validity in the conifold case 
of the holographic relation (12.631) . that can be derived only for the orbifold. In Appendix ICl 
we propose an alternative interpretation of the IR regime of our field theory, based on some 
nice observations made in [13] about the KW field theory. We argue that the theory may 
flow to a strongly coupled fixed point, although the string frame curvature invariant is large, 
as in the Klebanov-Witten solution for small values of g s N c . 

Contrary to the IR limit, the UV regime of the theory is dominated by flavors and we find 
the same kind of behavior for all values of the relevant deformations C\ and c<i- The gauge 
couplings increase with the energy, irrespective of the number of flavors. At a finite energy 
scale that we conventionally fixed to p = 0, the gauge theory develops a Landau pole, as told 
by the string coupling that diverges at that particular radius. This energy scale is finite, 
because p = is at finite proper distance from the bulk points p < 0. 

At the Landau pole radius the supergravity description breaks down for many reasons: 
the string coupling diverges as well as the curvature invariants (both in Einstein and string 

15 Here it is manifest why the SUGRA /3-function computed in this context with probe branes matches the 
field theory one, even if this requires the absence of order Nf /N c corrections to the anomalous dimensions 
7a,b, which one does not know how to derive (the stress-energy tensor is linear in Nf/N c ). It is because 
those corrections are really of order e^Nf/N c , and in the IR — > 0. 




09 




3N C - 2N C (1 - 7A ) - N f (l - 7/ ) 



(2.83) 
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frame), and the ip circle shrinks. An UV completion must exist, and finding it is an interesting 
problem. One could think about obtaining a new description in terms of supergravity plus 
branes through various dualities. In particular T-duality will map our solution to a system 
of NS5, D4 and D6-branes, which could then be uplifted to M-theory. Anyway, T-duality 
has to be applied with care because of the presence of D-branes on a non-trivial background, 
and we actually do not know how to T-dualize the Dirac-Born-Infeld action. We leave this 
interesting problem for the future. 

3 Part II: Generalizations 

In this section we are going to extend the smearing procedure of the D7-brane, which was 
formulated in section [2] for the particular case of the AdS$ x T 1 ' 1 space, to the more general 
case of a geometry of the type AdS§ x M5, where M5 is a five-dimensional compact manifold. 
Of course, the requirement of supersymmetry restricts greatly the form of M5. Actually, 
we will verify that, when M 5 is Sasaki-Einstein, the formalism of section [2] can be easily 
generalized. As a result of this generalization we will get a more intrinsic formulation of the 
smearing, which eventually could be further generalized to other types of flavor branes in 
different geometries. 

First of all we are going to generalize the effect of the smearing on the Wess-Zumino term 
of the D7-brane action for a general geometry. Following the line of thought that led to the 
action (11.61) and as is clear from the fact that it is linear in Cg, the smearing can be modelled 
by means of the substitution: 

Swz = T 7 J2 [ <% -> T 7 ! QAC 8 , (3.1) 

N J M s JMw 

where Q is a two-form which determines the distribution of the RR charge of the D7-brane 
in the smearing and A^io is the full ten-dimensional manifold. For a supersymmetric brane 
one expects the charge density to be equal to the mass density and, thus, the smearing of the 
DBI part of the D7-brane action should be also determined by the form Q. Let us explain in 
detail how this can be done. First of all, let us suppose that Q is decomposable, i.e. that it 
can be written as the wedge product of two one- forms. In that case, at an arbitrary point, Q 
would determine an eight-dimensional orthogonal hyperplane, which we are going to identify 
with the tangent space of the D7-brane worldvolume. A general two-form Q will not be 
decomposable. However, it can be written as a finite sum of the type: 

Q = Y" fi« , (3.2) 

where each fiM [ s decomposable. At an arbitrary point, each of the fiW's is dual to an eight- 
dimensional hyperplane. Thus, Q will determine locally a collection of eight-dimensional 
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hyperplanes. In the smearing procedure, to each decomposable component of Q we associate 
the volume form of its orthogonal complement in Aiio. Thus, the contribution of every 
to the DBI action will be proportional to the ten-dimensional volume element. Accordingly, 
let us perform the following substitution: 

Sdbi = -T 7 £ / d s £ yf^G~ s -> -T 7 [ d 10 x e* | fi« | , (3.3) 

N J Ms JMio 1 

where | | is the modulus of Vl^ and represents the mass density of the i th piece of Vl 
in the smearing. There is a natural definition of | fl® | which is invariant under coordinate 
transformations. Indeed, let us suppose that VL^ is given by: 

fiW = \\ S ^ N dx M Adx N . (3.4) 
Then, | | is defined as follows: 

| ^ W | = ^n® N n%G"rG»Q . (3.5) 

Notice that f2 acts as a magnetic source for the field strength F%. Actually, from the 
equation of motion of Cg one gets that Q is just the violation of the Bianchi identity for Fx, 
namely: 

dF-y = -Q . (3.6) 

For a supersymmetric configuration the form Q is not arbitrary. Indeed, eq. f !3.6p determines 
F\ which, in turn, enters the equation that determines the Killing spinors of the background. 
On the other hand, Q must come from the superposition (smearing) of K-symmetric branes. 
When the manifold M 5 is Sasaki-Einstein, we will show in Section 13.21 that Q can be deter- 
mined in terms of the Kahler form of the Kahler-Einstein base of M5 and that the resulting 
DBI+WZ action is a direct generalization of the result written in (11.61) . We will also show 
that the existence of Killing spinors implies that the functions appearing in the ansatz satisfy 
a system of first-order differential equations analogous to that written in (12.201) . 

3.1 General Smearing and DBI Action 

Here we will elaborate on the previous construction: writing the DBI action for a general 
smearing of supersymmetric D7-branes. We mean that in gen eral on an M = 1 background 
there is a continuous family of supersymmetric 4-manifolda 16 l that the D7-branes can wrap 

16 Even if we try to be general, we still stick to the case with vanishing B 2 background and vanishing Fmn 
on the brane world-volume. 
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corresponding to quarks with the same mass and quantum numbers. All these configurations 
preserve the same four supercharges, so we can think of putting D7's arbitrarily distributed 
(with arbitrary density functions) on these manifolds. We want to write the DBI plus WZ 
action for this system. 

Supersymmetry plays a key role. The fact that we can put D7's and not anti-D7's implies 
that the charge distribution completely specifies the system. For D7-branes the charge 
distribution is a 2-form Q, which can be localized (a "delta form" or current) or smooth (for 
smeared systems). The Bianchi identity reads gLF\ = — Q and is easily implemented through 
the WZ action ( 13. II) : Swz = T 7 j Q A Cg. Notice that a well defined Q not only must be 
closed (which is charge conservation) but also exact. Moreover the supersymmetry of this 
class of solutions forces Q to be a real (l,l)-form (with respect to the complex structure). 
Supersymmetry also guides us in writing the DBI action, because the energy distribution 
must be equal to the charge distribution. But there is a subtlety here, because the energy 
distribution is not a 2-form, and some more careful analysis is needed. 

Let us start considering the case of a single D7-brane localized on A4$. We can write its 
DBI action as a bulk lOd integral by using a localized distribution 2-form Q such that 

/ (fie* \[^&%= [ d l0 x e*>/=G|fi| • (3.7) 

Q is the Poincare dual to Ais- It can be (locally) written as Q = 5( 2 \M.g)\/ — G$/ V— GaA/5, 
through a properly normalized delta function and the product of two 1-forms (in general not 
separately globally defined) orthogonal to the 8-submanifoldJ^] In particular it is decompos- 
able. 

The decomposability of a 2-form can be established through Pliicker's relations, and the 
minimum number of decomposable pieces needed to write a general 2-form is half of its rank 
as a matrix^. So the decomposability of a 2-form at a point means that it is dual to one 8d 
hyperplane at that point; in general a 2-form is dual to a collection of 8d hyperplanes. 

If we do a parallel smearing of our D7-brane we get a smooth charge distribution 2- 
form, non-zero at every point. This corresponds to put a lot of parallel D7's and go to the 
continuum limit. Being the smearing parallel, we never have intersections of branes and the 
2-form is still decomposable. As a result (13 .7p is still valid. If instead we construct a smeared 
system with intersection of branes, the charge distribution Q is no longer decomposable. 
Every decomposable piece corresponds to one 8d hyperplane, tangent to one of the branes 
at the intersection. Since energy is additive, the DBI action is obtained by summing the 
moduli of the decomposable pieces (and not just taking the modulus of Q). Every brane 

17 This orthogonality does not need a metric. A 1-form is a linear function from the tangent space to R, 
and its kernel is a 9d hyperplane. The 8d hyperplane, tangent to the submanifold, orthogonal to the two 
1-forms, is the intersection of the two kernels. 

18 The rank of an antisymmetric matrix is always even. 
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at the intersection defines its 8d hyperplane and gives its separate contribution to the DBI 
action and to the stress-energy tensor. We simply sum the separate contributions because 
of supersymmetry: the D7's do not interact among themselves due to the cancellation of 
attractive/repulsive forces. Notice that in doing the smearing of bent branes, one generically 
obtains unavoidable self-intersections. 

Summarizing, given the splitting of the charge distribution 2-form into decomposable 
pieces f2 = J2k ^ > the DBI action reads 



S 



DBI 



The last step is to provide a well defined and coordinate invariant way of splitting the charge 
distribution Q in decomposable pieces. It turns out that the splitting in the minimal number 
of pieces compatible with supersymmetry is almost unique. 

In our setup, Q lives on the internal 6d manifold, which is complex and S , f/(3)-structure. 
This means that the internal geometry has an integrable complex structure X and a non- 
closed Kahler form J compatible with the metric: J^ — 9acJh- We can always find a 
vielbein basis that diagonalizes the metric and block-diagonalizes the Kahler form: 

^ (3.9) 

J = e 1 A e 2 + e 3 A e 4 + e 5 A e 6 . 

This pattern is invariant under the structure group SU(3) (without specifying the holo- 
morphic 3- form, it is invariant under U(3)), as is also clear by expressing them in local 
holomorphic basis: e Zl = e 2t ~ l + ie 2 \ e Zl = e 2i_1 — ie 2 \ with i = 1,2,3. One gets the 
canonical expressions: g = £\ e Zi ®s e Zi an d J = \c Zi A e z \ 

In our class of solutions, the supersymmetry equations force the charge distribution to 
be a real (1, l)-form with respect to the complex structure (see [12])- Notice that such a 
property is shared with J. The dilatino equation is e^Ff' 1 ^ = id<f) (which without sources 
amounts to the holomorphicity of the axio-dilation r = Co + i e~^). From this one gets 

n = ~dF 1 = 2i e-*(d(j) A 8(f)- dd(j)) . (3.10) 

It's manifest that Q is (1, 1) and Q* = Q. Going to complex components Q = QikC Zl A 
e Zk , the reality condition translates to the matrix being anti-hermitian. Thus it can 
be diagonalized with an SU(3) rotation of vielbein that leaves (13.91) untouched, and the 
eigenvalues are imaginary. 

Going back to real vielbein and summarizing, there is always a choice of basis which 
satisfies the diagonalizing condition ( 13. 91) and in which the charge distribution can be written 
as the sum of three real (1,1) decomposable pieces: 

Q = -Ai e 1 A e 2 — A 2 e 3 A e 4 — A 3 e 5 A e 6 . (3.11) 
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Supersymmetry forces the eigenvalues A a to be real and, as we will see, positive. Moreover, as 
inferred by the previous construction, the splitting is unique as long as the three eigenvalues 
A a are different, while there are embiguities for degenerate values, but different choices give 
the same DBI action. 

We conclude noticing that, in order to extract the eigenvalues \Xk\ = |fi( fc )| it is not 
necessary to construct the complex basis: one can simply compute the eigenvalues of the 
matrix (Q)mp9 PN in any coordinate basis. But in order to compute the stress-energy tensor, 
the explicit splitting into real (1,1) decomposable pieces is in general required. 



3.2 The BPS Equations for Any Sasaki-Einstein Space 



Let us now explain in detail the origin of the system of first-order differential equations (12.20p . 
As already explained in section [2], the system (I2.20p is a consequence of supersymmetry. 
Actually, it turns out that it can be derived in the more general situation that corresponds to 
having smeared D7-branes in a space of the type AdS§ x M5, where M5 is a five-dimensional 
Sasaki-Einstein (SE) manifold. Notice that the T 1,1 space considered up to now is a SE 
manifold. In general, a SE manifold can be represented as a one-dimensional bundle over a 
four- dimensional Kahler-Einstein (KE) space. Accordingly, we will write the M5 metric as 
follows 

ds 2 SE = ds 2 KE + (dr + A) 2 , (3.12) 

where d/dr is a Killing vector and ds 2 KE stands for the metric of the KE space with Kahler 
form J = dA / 2. In the case of the T 1,1 manifold the KE base is just S 2 x S 2 , where the 
S 2, s are parametrized by the angles (6*j, <fi) and the fiber r is parametrized by the angle if). 

Our ansatz for ten-dimensional metric in Einstein frame will correspond to a deformation 
of the standard AdS$ x M5. Apart from the ordinary warp factor h(r), we will introduce 
some squashing between the one form dual to the Killing vector and the KE base, namely: 



ds 2 



h(r) 2 dx\ 3 + h(r) 



dr 2 + 



^Us 2 KE + e 2f ^(dT + Af 



(3.13) 



Notice that, indeed, the ansatz (13.131) is of the same type as the one considered in eq. (11.91) for 
the deformation of AdS^ x T 1,1 . In addition our background must have a Ramond-Ramond 
five form: 

F 5 = K(r) dx° A ■ • • dx 4 A dr + Hodge dual, (3.14) 

and a Ramond-Ramond one-form Fi which violates Bianchi identity. Recall that this vio- 
lation, which we want to be compatible with supersymmetry, is a consequence of having a 
smeared D7-brane source in our system. Our proposal for F\ is the following: 



F 1 = C (dr + A) 



(3.15) 
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where C is a constant which should be related to the number of flavors. Moreover, the 
violation of the Bianchi identity is the following^: 

dF 1 = 2C J. (3.16) 

Notice that eq. (I3.16P corresponds to taking Q = —2CJ in our general expression (13.61) . To 
proceed with this proposal we should try to solve the Killing spinor equations by imposing 
the appropriate projections. Notice that the ansatz is compatible with the Kahler structure 
of the KE base and this is usually related to supersymmetry. 

Before going ahead, it may be useful for the interested reader to make contact with the 
explicit case studied in the previous section, namely the Klebanov-Witten model. In that 
case the KE base is 

ds KE = i£(d^ + sm 2 0,d V J) (3.17) 

i=l,2 

whereas the one form dual to the Killing vector d/dr is dr = dip/ 3 and the form A reads 

A = - ^ cos 6>! dip i + cos #2 d<£>2 j • (3.18) 

3 TV 

Moreover, the constant C was set to -r-*- in that case. 

' 4-7T 

Let us choose the following frame for the ten-dimensional metric: 

e*" = [h(r)]~* da? , e r = [ h{r) ] * dr , 



e 



°= [h(r)Y e f{r) (dr + A) e a = [ h(r) ] 1 e 9(r) e a 



(3.19) 



where e a a = 1, . . . , 4 is the one- form basis for the KE space such that ds 2 KE = e a e a . In 
the Klebanov-Witten model the basis taken in (12. 14j) corresponds to: 



e 1 = sin^x dipi , e 2 = dQ\ , 

e 3 = sin 9 2 dip 2 , e 4 = d9 2 ■ 

Let us write the five-form F 5 = JF 5 + *JF 5 of eq. (13.141) in frame components: 
T h = K{r) [h(r) } * e x ° A ■ • • A e x " A e r , 



(3.20) 



e 4 . 



*jr 5 = _ K ( r ) ^( r ) ] 4 go A . . . A g 4 = -Kh 2 e 49+f {dr + A) A e 1 A • • • A 

(3.21) 

The equation dF$ = immediately implies: 

Kh 2 e^ +f = constant = , (3.22) 

Vol(M 5 ) ' v ; 
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We are considering that J = ^J a bdx a Adx b and that the Ricci tensor of the KE space satisfies R a b = 6 g a b- 
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where the constant has been obtained by imposing the quantization condition (I2.16P for 
a generic M5. It will also be useful in what follows to write the one-form F\ in frame 
components: 

Fi = Ch-^e- f e°. (3.23) 
Let us list the non-zero components of the spin connection: 



Cj x r 



u 0r 



±h'h~U*\ (/i = 0,---,3), 

aj + 9']h^e\ (o=l,... ,4) 

r 1 h 



4l + // 



(3.24) 



- e^h-\j nh e b 



u" a = e- ~»n * j ab 
tu ab = u 



■ - ab - e f - 29 h-^ J ab e° 



where uj ab are components of the spin connection of the KE base. 

Let us now study under which conditions our ansatz preserves some amount of supersym- 
metry. To address this point we must look at the supersymmetric variations of the dilatino 
(A) and gravitino (ipu)- These variations have been collected in appendix [K\ for both the 
Einstein and string frame. We have written them in eq. (12.171) for the particular case in 
which the three-forms of supergravity are zero. Recall that the variations written in (12.171) 
correspond to the Einstein frame and we have used a complex spinor notation. 

It is quite obvious from the form of our ansatz for Fi in (13.231) that the equation resulting 
from the dilatino variation is: 

((f)' - ie^Ce' f r r0 ) e = . (3.25) 

In eq. (13.251) . and in what follows, the indices of the T-matrices refer to the vielbein compo- 
nents ( 13TT9|) . 

Let us move on to the more interesting case of the gravitino transformation. The space- 
time and the radial components of the equation do not depend on the structure of the internal 
space and always yield the following two equations: 

h' + K h 2 = , 

d r e - Ik he = . (3.26) 
8 

To get eq. (I3.26P we have imposed the D3-brane projection 

r^i^s e = — ie, (3.27) 
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and we have used the fact that the ten-dimensional spinor is chiral with chirality 

rx°...a: 3 r01234 € = € . (3.28) 

It is a simple task to integrate the second differential equation in (I3.26[) : 

e = h~h , (3.29) 

where e is a spinor which can only depend on the coordinates of the Sasaki-Einstein space. 

In order to study the variation of the SE components of the gravitino it is useful to 
write the covariant derivative along the SE directions in terms of the covariant derivative 
in the KE space. The covariant derivative, written as a one-form for those components, 
D = d + \ Cju r /J , is given by 

D = D — -J ab h-^e f - 29 T ab e° - - J ab /T* e f ~ 29 T ob e a + 
4 2 

+ \ h'\ (i h - + g>) r» e» + \ h ~ l {\j + f) T° r e° , (3.30) 

where D is the covariant derivative in the internal KE space. 

The equation for the SE components of the gravitino transformation is 

bit - Ik hi T rI e + -e+FjPe = 0. (3.31) 
8 4 

This equation can be split into a part coming from the coordinates in the KE space and a 
part coming from the coordinate which parameterizes the Killing vector. For this purpose, 
it is convenient to represent the frame one-forms e a and the fiber one-form A in a coordinate 
basis of the KE space 

e a = E a m dy m , 



(3.32) 
^4 = A m dy , 

with y m m — 1, . . . , 4 a set of space coordinates in the KE space. 

After a bit of algebra one can see that the equation obtained for the space coordinates y m 
is simply 

An e - \ J ab e 2(f - 9) A m T ab e-^ J ab h~$ e^ 9 E a m T ofe e + 
+ \h~- (\j + 9') E« m E- e + \ (\ \ + f) ef A m r 0r e - (3.33) 

- X -Kh* {E a m Y ra + h*e f A m T r0 )e + C A m e = , 
whereas the equation obtained for the fiber coordinate r is given by 



dr A JabC 6 + 2 [ A h + 1 )C 6 

- -Khe f T r0 e + -e^Ce = . 
8 4 



(3.34) 
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Let us now solve these equations for the spinor e. First of all, let us consider the dilatino 
equation (13. 25ft . Clearly, this equation implies that the spinor must be an eigenvector of the 
matrix T r o. Accordingly, let us require that e satisfies 

r r0 e — — ie . (3.35) 

Moreover, a glance at eqs. (I3.33f) and (13.341) reveals that e must also be an eigenvector of 
the matrix J a b^ ab - Actually, by combining eqs. (I3.27P , (13.281) and (13 .35|) one easily obtains 
that 

Ti 2 e = T 34 e . (3.36) 

To simplify matters, let us assume that we have chosen the one-form basis e a of the KE in 
such a way that the Kahler two-form J takes the canonical form: 

J = e 1 A e 2 + e 3 A e 4 . (3.37) 

In this basis, after using the condition ( 13.361) . one trivially gets: 

J ab T ab e = 4r 12 e . (3.38) 

Thus, in order to diagonalize J a b^ ab , let us impose the projection 

r 12 e = -ie , (3.39) 

which implies 

r 34 e=-ie, J ab T ab e = -Aie . (3.40) 

Let us now use the well-known fact that any KE space admits a covariantly constant 
spinor rj satisfying: 

3 

D mV = —ziA m r] , (3.41) 
from which one can get a Killing spinor of the five-dimensional SE space as: 

e = e~ l ^ T T] . (3.42) 



Actually, in the KE frame basis we are using, rj turns out to be a constant spinor which 
satisfies the conditions T 12 rj = r 34 r] = —it]. Let us now insert the SE Killing spinor e of 
eq. ( I3.42p in our ansatz (13.291) . i.e. we take the solution of our SUSY equations to be: 



3 , 



e = h~s e-2 lT r] . (3.43) 

By plugging (I3.43j) into eqs. ( 13.331) and (I3.34[) . and using the projections imposed to e and 
(I3.4ip . one can easily see that eqs. (13.331) and (13.341) reduce to the following two differential 
equations: 

l^ + g'+lKh-e^ = 0, 
Ah 4 

~ + /' + \Kh + 2e f - 29 - 3e~ f + -e+~ f = . (3.44) 
4 h 4 2 
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By combining all equations obtained so far in this section we arrive at a system of first-order 
BPS equations for the deformation of any space of the form AdS^ x M5: 

6' - Ce*- f = , 



9 - e 



Vol(M 5 ) 
f-2g = o 



. 



/' +2e f ~ 29 - 3e~ f + -e^ f = . (3.45) 

Notice that, indeed, this system reduces to the one written in eq. (12.201) for the conifold, if we 
take into account that for this later case the constant C is 3Nf/(4ir) and Vol{T 1,1 ) = 167r 3 /27. 

It is now a simple task to count the supersymmetries of the type (I3.43P preserved by our 
background: it is just thirty-two divided by the number of independent algebraic projection 
imposed to the constant spinor rj. As a set of independent projections one can take the ones 
written in eqs. (I3.27p . (I3.35P and (13.391) . It follows that our deformed background preserves 
four supersymmetries generated by Killing spinors of the type displayed in eq. (13.431) . 



3.3 The BPS and Einstein Equations 

In this section we will prove that the BPS system implies the fulfilment of the second-order 
Euler-Lagrange equations of motion for the combined gravity plus brane system. To begin 
with, let us consider the equation of motion of the dilaton, which can be written as: 

d M (G MN V^G d N <p) = Fl-^^- S DBI , (3.46) 



where Gmn is the ten-dimensional metric. Using the DBI action (13. 3p for the smeared 
D7-branes configuration, we find: 

-^£*»=*£,l° w l- < 3 ' 47 > 

The charge density distribution is Q = —2CJ (see eq. (I3.16P ). Recall that the Kahler form 
J of the KE base manifold has the canonical expression (13.371) . It follows that Q has two 
decomposable components given by: 

= -2C e 1 A e 2 = -2C e~ 29 e 1 A e 2 , 

1 (3-48) 
= -2C e 3 A e 4 = -2C e ~ 29 e 3 A e 4 , 
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where the e a one-forms have been defined in (I3.19p . Therefore, the moduli of the fiW's can 
be straightforwardly computed: 

= |fi( 2 )| = 2\C\h~^e- 29 . (3.49) 

By using the explicit form of the metric, our ansatz for F\ and the previous formulae (13.491) 
one can convert eq. (13.461) into the following: 

0" + (V + f) 0' = C 2 e 2 ^ 2/ + 4 |C| e^ 29 . (3.50) 

It is now a simple exercise to verify that eq. (I3.50P holds if the functions <p, g and / solve 
the first-order BPS system (I3.45P and the constant C is non-negative. In what follows we 
shall assume that C > 0. 

To check the Einstein equation we need to calculate the Ricci tensor. In flat coordinates 
the components of the Ricci tensor can be computed by using the spin connection. The 
expression of the curvature two-form in terms of the spin connection is 

Rmn = d&MN + ^MP^^N ' (3.51) 

with the curvature two-form defined as follows: 

By using the values of the different components of the ten-dimensional spin connection 
written in (I3.24p we can easily obtain the Riemann tensor and, by simple contraction of 
indices, we arrive at the following flat components of the Ricci tensor: 



1 



2 



1 



lh" l(h'\ lh', h' 



4(</) 2 - g'f - 2e 2/ ~ 49 + 6e" 29 



R m = , M^N . 



(3.53) 
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From these values it is straightforward to find the expression of the scalar curvature, which 
is simply 



, i / lh" lh' „, h' „ , /N9 

R = -h-^-- + ~ 1 -f> + 2 1 - 9 > + 8 g » + 20( g y + 

+ 8( ? 7' + 2f / + 2(f) 2 + 4e 2/ - 4s -24e- 2 ^ . (3.54) 

Let us evaluate the different contributions to the right-hand side of Einstein's equations. 
The contributions from the five- and one-forms have been written in the first equation in 
(11. 7p and is immediately computable from our ansatz of eqs. (13. 14j) and A3. 15j) . On the other 
hand, the contribution of the DBI part of the action is just 

2k\ 5S DB i , s 

- 7^ Jg^ ■ ^ 

By using our expression ( 13. 3ft of Sdbi, with Q = —dFi, together with the definition (13.51) . 
one easily arrives at the following expression of the stress-energy tensor of the D7-brane: 



T- - = -— 

± MN 2 



where we have used that 2k1 Tj = 1 and we have written the result in flat components. By 
using in (13.561) the values given in eqs. 0338) and Q325D of and its modulus, we arrive 
at the simple result: 

^b*^ = ~2C h 2 s ^ x i x j , 
T rr = Too = -2Ch~^e*- 29 , 



T ab = -C /T5 e^ 29 5 ab , (a, 6 = 1, • • ■ , 4) , (3.57) 

where the indices refer to our vielbein basis ( 13. 19ft . As a check of this result one can explicitly 
verify that the result of eq. (I2.2ip for the conifold, when written in flat indices, reduces to 
the simple expressions written in (13.57p . 

With all this information we can write, component by component, the set of second order 
differential equations for h, g, f and <ft that are equivalent to the Einstein equations. One can 
then verify, after some calculation, that these equations are satisfied if <fi and the functions of 
our ansatz solve the first-order system ( 13.45ft . Therefore, we have succeeded in proving that 
the background obtained from the supersymmetry analysis is a solution of the equations 
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of motion of the supergravity plus Born-Infeld system. Notice that the SUSY analysis 
determines Fx, i.e. the RR charge distribution of the smeared D7-branes. What we have just 
proved is that eq. (13.561) gives the correct energy-momentum distribution associated to the 
charge distribution Q = —dFx of smeared flavor branes. 

To finish this section let us write the DBI action in a different, and very suggestive, fashion. 
It turns out that, for our ansatz, the on-shell DBI action can be written as the integral of a 
ten-form and the corresponding expression is very similar to the one for the WZ term (eq. 
(13. ip ). Actually, we show below that 

Sdbi = -T 7 ! e^QAQs, (3.58) 

JMio 

where f2§ is an eight-form which, after performing the wedge product with the smearing 
two-form Q, gives rise to a volume form of the ten-dimensional space. Let us factorize in f2§ 
the factors coming from the Minkowski directions: 

tt s = ft-'^A^, (3.59) 

where Q4 is a four-form in the internal space. Actually, one can check that Q4 can be written 
as: 

n 4 = ^JAJ, (3.60) 

where J is the following two-form: 

J = hSe 29 J + hh f dr A (dr + A) . (3.61) 

To verify this fact, let us recall that Q = 2CJ and thus 

n Afi 8 = Ch- l d A x A J A J A J . (3.62) 

Taking into account that | J A J is the volume form of the KE base of M 5 , we readily get: 

d 4 x A J A J A J = Ae- 29 h* V^G d w x , (3.63) 

from where one can easily prove that eq. (I3.58P gives the same result as in equation (13.31) 
with Q = —dFx. 

3.4 A Superpotential and the BPS Equations 

It is interesting to obtain the system of first-order BPS equations (I3.45P by using an alterna- 
tive approach, namely by deriving them from a superpotential. Generically, let us consider 
a one-dimensional classical mechanics system in which rj is the "time" variable and A(rj), 
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$ m (?7) (m = 1, 2 . . .) are the generalized coordinates. Let us assume that the Lagrangian of 
this system takes the form: 



-A 



k (d v A) 2 - -G mn ^)d^ m d^ n - V{$) 



(3.64) 



where k is a constant and V($) is some potential, which we assume that is independent of 
the coordinate A. If one can find a superpotential W such that: 



y ' 2 <9$ m d$" 4k 



(3.65) 



then the equations of motion are automatically satisfied by the solutions of the first order 
system: 

A J i A<\>>" nw 

(3.66) 



dA 

drj 



2k 



dW 



dr) <9$ n 

Let us now show how we can recover our system (I3.45P from this formalism. The first step 
is to look for an effective Lagrangian for the dilaton and the functions of our ansatz whose 
equations of motion are the same as those obtained from the Einstein and dilaton equations 
of Type IIB supergravity. One can see that this lagrangian is: 

1 



L 



eff 



R 



h-* (0') 



'\2 



(3.67) 



(3.68) 



where R is the scalar curvature as written in (I3.54j) and Q is the constant 

{2^N C 
W Vol(M 5 ) ' 

The Ricci scalar (13.54p contains second derivatives. Up to total derivatives L e ff takes the 
form: 

2 



J eff 



-Ag+f 



1 / ti 



t II + 12 (g') 2 + 8g'f - Ae 2f - ig + 24 e~ 2a 



Ll\2 



2 \ h 



-91 u-K 

2 



-2/ 



c 1 



e 2(<p-f) _ ACe 4>-2 g 



(3.69) 



We want to pass from the lagrangian (I3.69|) to that in eq. (I3.64p . With that purpose in mind 
let us perform the following redefinition of fields: 



AA 



h*e 



:P-9 



hSe 2 \ 



3 2/ 



hh 2f . 



In addition, we need to do the following change of the radial variable lj 

dr 
drj 



A_8 

e 4 3 



9-1/ 



(3.70) 



(3.71) 



20 The change of the Lagrangian under that change of radial variable is L e ff = ^ L e ff. 
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Once we have done the previous redefinitions, the Lagrangian we obtain is: 



J eff 



(3.72) 



where the dot means derivative with respect to 77 and V(g, f, 0) is the following potential: 

\-(<i.f.o) = <-i<&+f)[ 4e 2/ ~ 4 ^ -2Ae~ 2 ' 9 + 5! e -m+f) + 91 + AC 

2 2 

(3.73) 

The above lagrangian has the desired form (see eq. ( 13.641) ) and we can identify the constant 
k and the elements of the kinetic matrix G m „ as: 



3 
4 



56 

'->' = - • Ggg = y , Gff = 3 , Ggf = 77 > G 1 



8 
3 



(3.74) 



One can now check that, given the above expression of the potential, the following superpo- 
tential 



W 



4(43+/) 



Q e -45-f - 4e f - 2 ' 9 - Qe- f + Ce< 



-f 



(3.75) 



satisfies eq. ( I3.65P for the values of k and G mn written in eq. (13.741) . It is now immediate to 
write the first-order differential equations that stem from this superpotential. Explicitly we 
obtain: 



A 



W 



f 



-e 3 



1(49+/) 



4 6 3 



1(49+/) 



Q e -*9-f + 4e /-2§ 

-Qe-^-f - 8e f '- 2 ' 9 + 12 e~ f - 2Ce^~ f 



C e <t>-i(9+f) _ 



(3.76) 



In order to verify that this system is equivalent to the one obtained from supersymmetry, let 
us write down explicitly these equations in terms of the old radial variable (see eq. (13.711) ) 
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and fields (see eqs. (13.701) ). One gets: 
h' 

— + 8 a' + 2f = -Qhr l e- Ag -f + 4e / ~ 29 + 6e _/ - Ce*~ f , 
h 

+ = t f - 29 -\Qh-ie-*-* , 

™ + /' = 3e~f - 2e'" 2 * - \q h~ x e~ i9 ~ f - -C , 
Ah 4 2 

0' = Ce^ , (3.77) 

which are nothing else than a combination of the system of BPS equations written in (I3.45p . 

Let us now use the previous results to study the 5d effective action resulting from the 
compactification along M5 of our solution. The fields in this effective action are the func- 
tions / and g, which parameterize the deformations along the fiber and the KE base of M 5 
respectively, and the dilaton. Actually, in terms of the new radial variable r\ introduced in 
(13.7ip . the ten-dimensional metric can be written as: 

ds 2 = e"i (/+4§) leidx^dx^ + d V 2 } + e 2 ' 9 ds 2 KE + e 2f ~ (dr + A) 2 . (3.78) 



The corresponding analysis for the unflavored theory was performed in [81 [35] . For simplicity, 
let us work in units in which the AdS^ radius L is one. Notice that the quantity Q defined 
in (I3.68P is just Q = 4L 4 . Thus, in these units Q = 4. To make contact with the analysis 
of refs. El [35] , let us introduce new fields q and p which, in terms of / and g are defined as 
followS 

q=l(f + 4g), P=-\U-g). (3.79) 

15 5 

In terms of these new fields, the potential (13.73P turns out to be 

V(p,q,<f>) = Ae- Sq ~ 12p - 24e- 8q - 2p + — e 2 *~ 8q+8p + 8e~ 20q + 4Ce (f> - 8q ~ 2p , (3.80) 
and the effective lagrangian f)3.72p can be written as: 

L eff = V=g~s[R 5 - 1<P 2 - 20p 2 - 30q 2 - V] , (3.81) 
where g§ = —e 2A is the determinant of the five-dimensional metric ds 2 , = dx^dx^ + drf 



and i? 5 = -\2A+\A 2 



is its Ricci scalar. One can check that the minimum of the potential 



(I3.80p occurs only at p = q = = 0, which corresponds to the conformal AdS§ x M 5 



21 



The function p is called / in refs. [8\ [35] . 
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geometry. Moreover, by expanding V around this minimum at second order we find out that 
the fields p and q defined in (13.791) diagonalize the quadratic potential. The corresponding 
masses are m 2 = 12 and m 2 = 32. By using these values in the mass-dimension relation 
( l2~80l) . we get: 

ml = 12 A„ = 6 

2 . (3.82) 

m 2 = 32 =>• A 9 = 8 . 

These scalar modes p and q are dual to the dimension 6 and 8 operators discussed in section 

El 



3.5 General Deformation of the Klebanov-Witten Background 

In this section we will explore the possibility of having a more general flavor deformation of 
the AdSs x T 1,1 background. Notice that, as T ' is a U(l) bundle over S 2 x S 2 , there exists 
the possibility of squashing with different functions each of the two S 2, s of the KE base. In 
the unflavored case this is precisely the type of deformation that occurs when the singular 
conifold is substituted by its small resolution. For this reason, it is worth considering this 
type of metric also in our flavored background. To be precise, let us adopt the following 
ansatz for the metric, five-form and one-form: 

ds 2 = h~ 1/2 dxl 3 + h 1/2 \dr 2 + ^(d6 2 + sin 2 0< dip 2 ) + —(d^ + cos ^ d ¥i) 2 \ 

^ i=l,2 ' i=l,2 * 



'1 + *)d 4 xAKdr 



C 

F\ = —(dif) + cos 9 2 dip 2 + cos Oxdifx) , 

(3.83) 

where C = 3iVy/47r, all functions depend on r and <7i(r) and #2( r ) are, in general, different 
(if gx = #2 — 9 we recover our ansatz (11.91) ). The equation dF 5 = immediately implies: 

Kh 2 e 2g 1+ 2g 2 +f = 2 ^N C = Q , (3.84) 

which allows to eliminate the function K in favor of the other functions of the ansatz. 
By following the same steps as in the gx = gi case and requiring that the background 
preserve four supersymmetries, we get a system of first-order BPS equations for this kind of 
deformation, namely: 

(f)' = Ce+~ f , 

h' = -Qe'f' 29 ^ 2 ^ , 

S&=e?-*», (< = 1,2), ( 3 - 85 ) 

f = 3 e -f - e f - 2 ^ - e f ~ 292 - - e^ f . 
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Notice that, as it should, the system (I3.85j) reduces to eq. (I3.45j) when g t = g 2 . 

It is not difficult to integrate this system of differential equations by following the same 
method that was employed for the g\ = g 2 case. First of all, we change the radial coordinate: 

dr = e f dp , (3.86) 

what allows us to get a new system: 

<P= Ce* , 

h= -Q e - 291 ~ 292 , 

^=e 2 ^, (. = 1,2), ( 3 - 87 ) 

/ = 3 - e 2/ " 291 - e 2f ~ 292 -je*, 

where now the derivatives are taken with respect to the new variable p. 

The equation for the dilaton in ( 13. 87ft can be integrated immediately, with the result: 

= -i - , (p < 0) , (3.88) 

u p 

where we have absorbed an integration constant in a shift of the radial coordinate. Moreover, 
by combining the equations for gi and g 2 one easily realizes that the combination e 2gi — e 292 
is constant. Let us write: 

On the other hand, by using the solution for <j)(r) just found and the equations for the g^'s 
in ( 13.871) . the first-order equation for / can be rewritten as: 

/ = 3 - gi - g 2 + ±- , (3.90) 

ip 

which can be integrated immediately, to give: 

with c being an integration constant. This constant can be absorbed by performing a suitable 
redefinition. In order to make contact with the case in which g\ = g 2 let us take c = 6. Then, 
by combining ( 13.911) with the equation of g 2 , we get 

e 4 92 +2 9l = e 2 gi+ 2 92+ 2f = _ Qf)e 6p ^ ( 3>92 ) 

which, after using the relation ( 13.891) . can be integrated with the result 

e 692 + ^a 2 e 492 = (1 - 6p) e 6p + a . (3.93) 
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Notice that, indeed, for a = this equation reduces to the gi = g 2 solution (see eq. (I2.37P ). 
Moreover, by combining eqs. ( 13.891) and ( 13.911) the expression of / can be straightforwardly 
written in terms of g 2 , as follows: 



,2/ 



(3.94) 



+ c 2 . 



(3.95) 



e 4 3 2 + a 2 e 2 ^ 

It is also easy to get the expression of the warp factor h: 

Kp) = ~ Q J e^+a 2 e 2 ^ 

Thus, the full solution is determined in terms of e 2g2 which, in turn, can be obtained from 
( I3.93P by solving a cubic algebraic equation. In order to write the explicit value of e 2g2 , let 
us define the function: 

Then, one has: 



= (1 - 6p) e 6 " + d . 



^92 



1 

2 



a 2 + 



where the function ((p) is defined in terms of £(p) as: 



r + [C(p)]' 



(3.96) 
(3.97) 



C(P) = 4C(p) 



4\/^(p) 2 -y£(p)- 



(3.98) 



In expanding these functions in series near the UV (p — > 0) one gets a similar behavior to 
the one discussed in section 12.31 Very interestingly, in the IR of the field theory, that is when 
p — > — oo, we get a behavior that is "softened" respect to what we found in section [2731 This 
is not unexpected, given the deformation parameter a. Nevertheless, the solutions are still 
singular. Indeed, the dilaton was not affected by the deformation a. 



3.6 Massive Flavors 

In the ansatz we have been using up to now we have assumed that the density of RR charge 
of the D7-branes is independent of the holographic coordinate. This is, of course, what is 
expected for a flavor brane configuration which corresponds to massless quarks. On the 
contrary, in the massive quark supersymmetric D7-brane has a non-trivial profile 

in the radial direction [20] and, in particular ends at some non-zero value of the radial 
coordinate. These massive embeddings have free parameters which could be used to smear 
the D7-branes. It is natural to think that the corresponding charge and mass distribution 
of the smeared flavor branes will depend on the radial coordinate in a non-trivial way. 

It turns out that there is a simple modification of our ansatz for Fi which gives rise to a 
charge and mass distribution with the characteristics required to represent smeared flavor 
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branes with massive quarks. Indeed, let us simply substitute in (13.451) the constant C by a 
function C(r). In this case: 

Fi = C{r) (dr + A) , 
dF 1 = 2 C{r) J + C'{r)dr A (dr + A) . y ' J 

Notice that the SUSY analysis of sect. 13.21 remains unchanged since only F\, and not its 
derivative, appears in the supersymmetric variations of the dilatino and gravitino. The final 
result is just the same system (13.451) of first-order BPS equations, where now one has to 
understand that C is a prescribed function of r, which encodes the non-trivial profile of the 
D7-brane. Notice that C(r) determines the running of the dilaton which, in turn, affects the 
other functions of the ansatz. 

A natural question to address here is whether or not the solutions of the modified BPS 
system solve the equations of motion of the supergravity plus branes system. In order to 
check this fact, let us write the DBI term of the action, following our prescription (13. 3p . 
Notice that, in the present case, Q = —dF\ is the sum of three decomposable pieces: 

Q = tt w + Q {2) + fi (3) , (3.100) 

where f^ 1 -* and f^ 2 -* are just the same as in eq. (13.481) . while f^ 3 -* is given by: 

fi( 3 ) = —C'(r) dr A (dr + A) = -h~* e' s C'(r) e r A e° . (3.101) 

The modulus of this new piece of Q can be straightforwardly computed, namely: 

|J) (3) | = /T3 e~ f I C'{r) | . (3.102) 

By using this result, together with the one in (13.491) . one readily gets the expression of the 
DBI terms of the action of the smeared D7-branes: 



Sdbi — —T-j / h 



C(r) | e~ 29 + | C\r) \ e~ f ) V^G d 10 x . (3.103) 



'Mio 

From this action it is immediate to find the equation of motion of the dilaton, i.e.: 

<J>" + (V + /') 0' = C 2 e 2 *" 2/ + 4 |C| e*~ 29 + \C'\. (3.104) 

It can be verified that the first-oder BPS equations (13.451) imply the fulfilment of eq. (13. 1Q4[) . 
provided the functions C(r) and C'(r) are non-negative. Notice that now, when computing 
the second derivative of from the BPS system (13.451) with C = C(r), a new term containing 
C'(r) is generated. It is easy to verify that this new term matches precisely the last term on 
the right-hand side of (13.1041) . 

It remains to verify the fulfilment of Einstein's equation. The stress-energy tensor of the 



brane can be computed from eq. (13.561) . where now the extra decomposable piece of dF\ 
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must be taken into account. The result one arrives at, in the vielbein basis (I3.19p . is a direct 
generalization of (13.571) : 



- x 1 x^ 



-e h~* 



T, 



ab 



2\C(r)\e- 29 + -\C'(r)\e- f \ W , (i, j = 0, . . . , 3) 



C(r) | e~ 29 + i | C"(r) | e" / 



(a, 6= 1,...,4) 



T 



00 



2 | C(r) | /T* e 4 



-2g 



(3.105) 



As happened for the equation of motion of the dilaton, one can verify that the extra pieces 
on the right-hand side of (I3.105P match precisely those generated by the second derivatives 
appearing in the expression (13.531) of the Ricci tensor if C(r) and C'(r) are non-negative. As 
a consequence, the first-order equations (13.451) with a function C{r) also imply the equations 
of motion for the ten-dimensional metric gMN- It is also interesting to point out that, if C(r) 
and C'(r) are non-negative, Sdbi can also be written in the form (13.581) . where Qs is exactly 
the same eight-form as in eqs. (13 . 59[) and (13.601) . 

Notice that, if the function C(r) = 3Nf(r)/4ir has a Heaviside-like shape "starting" at 
some finite value of the radial coordinate, then our BPS equations and solutions will be the 
ones given in section 12.31 for values of the radial coordinate bigger than the "mass of the 
flavor" . However, below that radial value the solution will be the one of Klebanov-Witten 
(or deformations of it, see appendix B), with a non-running dilaton. Aside from decoupling 
in the field theory, this is clearly indicating that the addition of massive flavors "resolves" 
the singularity. Physically this behavior is expected and makes these massive flavor more 
interesting. 



4 Summary, Future Prospects and Further Discussion 

In this paper we followed the method of [TF] to construct a dual to the field theory defined by 
the Klebanov and Witten after Nf flavors of quarks and antiquarks have been added to both 
gauge groups. In Section [2] of this work, we wrote BPS equations describing the dynamics 
of this system and found solutions to this first order system, that of course solves also all 
the second order equations of motion. We analyzed the solutions to the BPS system and 
learnt that, even when singular, the character of the singularity permits to get field theory 
conclusions from the supergravity perspective. 

We proposed a formulation for the dual field theory to these solutions, constructing a 
precise 4-dimensional superpotential. We studied these solutions making many matchings 
with field theory expectations that included the R-symmetry breaking and Wilsonian beta 
function. Also, using the well known (supergravity) superpotential approach, we learnt that 
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our field theory, aside from being deformed by a marginal (then turned irrelevant) operator, 
modifies its dynamics by giving VEV to operators of dimension six and eight. We explained 
how to change relations between couplings and ^-angles in the theory, from the perspective 
of our solutions. We believe that these many checks should encourage other physicists to 
study this background more closely. 

In Section [3] of this paper, we presented a careful account of the many technical details 
regarding the derivation of the results in section [2] summarized above. But most interestingly, 
Section [3] is not only about technical details. Indeed, using the logic and intuitions developed 
in Section [21 we generalized the approach described there for any five dimensional manifold 
that can be written as a Sasaki-Einstein space (a one- dimensional fibration over a Kahler- 
Einstein space). It is surprising that the same structure of BPS eqs and ten-dimensional 
superpotential repeats for all the manifolds described above. This clearly points to some 
"universality" of the behavior of 4-dimensional M = 1 SCFT's with flavors. 

We have added some brief comments about what happens when we take the number of 
flavors Nf = in our BPS eqs (see Appendix [Bj . It is interesting to recover some solutions 
studied in the past from this perspective since it puts into context previous analysis. Again, 
the careful study of this "unflavored" solutions might be of interest to many physicists. We 
shortly commented on the possibility of adding to the dynamics of the 4-d field theory fun- 
damentals with mass, presenting a general context to do this. We will exploit this procedure 
in the future to get a better understanding of our singular backgrounds, make contact with 
field theory results and study many other interesting problems. 

All the results described above not only motivate a more detailed analysis of this approach 
from a field theoretical viewpoint, but also emphasize the need for a deeper geometrical study, 
that clearly will reveal interesting underlying structure. 

4.1 Future Directions 

Many things can be done following the results of this paper. It is natural to extend the 
treatment to the case of the Klebanov-Tseytlin and Klebanov-Strassler solutions. The result 
is likely to be interesting, since the fundamentals and the KT cascade "push in different 
directions" in the RG flow. One might find a fine-tuned situation in which the IR dynamics 
is different from the one in the Klebanov-Strassler model. 

Other things that immediately come to mind are to study the dynamics of moving strings 
in this backgrounds, details related to dibaryons, flavor symmetry breaking, etc. Even when 
technically involved, it should be nice to understand the backreaction of probes where the 
worldvolume fields have been turned on, since some interesting problems may be addressed. 

The formalism developed to deal with configurations of IIB dual to massive fundamentals 
seems useful in different contexts. Needless to say, the approach adopted here is immediately 
generalizable to the case of type IIA backgrounds. Duals to M = 1 field theories have been 
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constructed and it seems natural to apply our methods in those cases. 

Finding black hole solutions in our geometries is not an elementary task; but it should 
not be very difficult. The interest of this problem resides in the fact that this will produce 
a "well-defined" black hole background where to study, among other things, plasmas that 
include the dynamics of color and flavor at strong coupling. This is a very well defined 
problem that we believe of much interest. 

On the field theory side, it should be interesting to understand in more detail how the 
smearing procedure affects the superpotential. We gave a possible answer and detailed 
study can uncover interesting subtleties. Here again, similar ideas can be extended to other 
situations in type IIA and type IIB. Getting a better handle on the field theory interpretation 
of our "generalized" approach of Part II seems also interesting. Indeed, understanding in 
detail what is the "universality" that produces the same dynamics for a large class of M = 1 
SCFT's with flavor would be nice. 

4.2 Further Discussion 

Let us finish this paper with some discussions that might be of interest for the reader. The 
first point we want to address is what could be the application of these results to Physics. 
Indeed, it is not easy to find an interesting physical system displaying a Landau pole (without 
a UV completion, like QED has, for example). Of course, as explained above, this paper 
is a first step in a more detailed study of a cascading field theory with flavors, that with 
no doubt has applications in Physics. Nevertheless, one can find some interesting problems 
already at this stage. 

As described above, finding a black hole in our geometry, might be a good simulation of 
the Physics of a strongly coupled quark-gluon plasma. Even more, since we would be only 
interested in effects in the hydrodynamics regime, using the IR of this black hole solutions 
should be enough to learn about Physics at RHIC, for example. One can also think that 
our paper starts the study of the different phases of this generalized Af = 1 SQCD obtained 
from Klebanov-Witten-like models. 

Let us change the subject of the discussion and go back to our procedure, that was well 
explained in the introduction of this paper. The reader may remember the difference between 
a weakly gauged symmetry and a global symmetry, let us now connect this to supergravity. 
One important distinction between the approaches for finding string duals to field theories 
with fundamentals is that in the approach where the solution consists only of supergravity 
fields, the field theory will have this global symmetry weakly gauged. On the contrary, in our 
case backreacting with the Born-Infeld action, the symmetry will be global. We can see this 
clearly in the fact that the BI action has the freedom to add worldvolume gauge fields (and 
scalars), hence introducing a gauged symmetry in the bulk, dual to a global symmetry in 
the boundary. In a (complicated) reduction of our Type IIB plus Born-Infeld action to five 
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dimensions, we would see some SU(Nf) gauge fields (as many as branches of flavor branes 
we added) that would enter in the holographic formulas to compute field theory correlators. 

It is interesting to notice that depending on the physical situation we want to work with, 
we should choose the approach used here or the complementary one of finding a solution 
purely in supergravity. Indeed, for situations where we do not want to take into account 
the "flavor degrees of freedom" of the extra branes, but what we want is to introduce some 
operator in the dual field theory (like a giant graviton, a domain wall or a Wilson line) we 
should work within the purely supergravity approach [43] . Indeed, if we are thinking about 
the presence of an operator (say in N = 4 SYM), there should be no "flavor degrees of 
freedom" in the solutions. 

Finally, we would like to comment on the smearing procedure. One way in which we can 
think about it is to realize that usually (unless they are D9 branes) the "localized" flavor 
branes will break part of the isometries of the original background dual to the unflavored field 
theory. The "smeared" flavor branes on the other hand reinstate these isometries (global 
symmetries of the field theory dual). In some sense the flavor branes are 'deconstructing' 
these dimensions (or these global groups) for the field theory of interest. In the case in which 
we have a finite number of flavors, these manifolds become fuzzy, while for Nf — ► oo, we 
recover the full invariance. 
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A Appendix: SUSY Transformations in String and 
Einstein frame 

The supersymmetry transformations of Type IIB supergravity were found long ago in ref. 
[44] . Here we will follow the conventions of the appendix A of [45], where they are written 
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in string frame. Let us recall them: 



1 1 

/-;;;; v/ ,].— ^ jr- 

(a.i; 



_i_ i±p( 3 ) p(s) MNP \ (s) 
2^] r MNP L °l) e 



+ ^4V (S)7VPQ -! + ^i&«rr«^ OT (<a a ))r« c w , 

where the superscript s refers to the string frame, <7j % = 1, 2, 3 are the Pauli matrices, if 
is the NSNS three-form and and F^ 5 ) are the RR field strengths. In (IA.10 e is a 

doublet of Majorana-Weyl spinors of positive chirality. 

We can study how these equations change under a rescaling of the metric like 

9mn = e * 9mn- (A.2) 

In doing that it is useful to follow Section 2 of [IB]. Under the above change for the metric, 
there are some quantities which also change: 

ry = e * T M , 
e« = ef e , 
A« = e"s A , 

^ = e-f(^-irWA«). (A.3) 

The equation for the dilatino in the new frame can be easily obtained whereas in doing the 
same for the gravitino equation we will use that 



V A /e + ^r A ?(Vjv0) + l(y M <f>) . (A.4) 



After some algebra with gamma-matrices, the SUSY transformations in Einstein frame we 
obtain are the following ones: 

+^4 5 Vr w (- 2 )r A/e . (A.5) 

In order to write the expression of the SUSY transformations, it is convenient to change 
the notation used for the spinor. Up to now we have considered the double spinor notation, 
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namely the two Majorana-Weyl spinors e\ and e 2 form a two-dimensional vector ( . We 

n ^ 

can rewrite the double spinor in complex notation aa 22 l e = 61 — 262. It is then straightforward 
to find the following rules to pass from complex to real spinors: 

e* <-> cr 3 e , —ie* <->• a - ! e , ie <-> za 2 e , (A. 6) 

where the Pauli matrices act on the two-dimensional vector 




B Appendix: The Unflavored Solutions 

In this appendix we will study our BPS system of linear ordinary differential equations (12.301) - 
( 12.33!) an d we will find its general solution in the absence of D7-branes. Not only we will 
recover the solution describing a stack of D3-branes placed at the apex of the real Calabi-Yau 
cone over a generic Sasaki-Einstein 5-manifold M 5 , preserving (at least) four supercharges, 
and its near-horizon limit AdS^ x M 5 dual to the (at least) Af = 1 superconformal gauge 
theory describing the IR dynamics on the stack of D3-branes, but also the solution describing 
D3-branes smeared homogeneously on a blown-up 4-cycle inside the Calabi-Yau, discussed 
in the paper [31] for the case of the conifold (more precisely a Z 2 orbifold of it) and then in 
full generality in [35] for all Calabi-Yau cones. We will also study the unflavored limit of the 
general deformation of the KW model analyzed in section 13.51 and we will show that it gives 
rise to the two-parameter metrics found in ref. [31] . 

Let us look at our BPS system of linear ordinary differential equations (I2.30H2.33]) . We will 
sometimes refer to the case of the conifold for the sake of simplicity. The generalization to 
any Sasaki-Einstein is straightforward, the only difference being the normalization in (12.331) 
and in the RR 5-form field strength, related to the volume of the Sasaki-Einstein base. 

First of all notice that Nf must be set to zero in the system of equations and not in our 
solution, since when we solved the equation for the dilaton (I2.32j) we supposed that Nf ^ 0. 
This allowed us to get ( 12.341) after shifting the radial variable. 

It is easy to show that the most general solution to the BPS system when Nf = 0, up to 

22 Notice that there is an ambiguity in the choice of the relation between complex and real spinors. 
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redefinition of the coordinates, is the following: 

00) = 0o (B.l) 

,9(p)= [e 6 " + Cl ] 1/6 (B.2) 

/(p) = e 3 ^e 6 " + Cl ]- 1/3 (B.3) 



e 



h(p) = c 2 - 4L 4 J dp e~ 49{p) (B.4) 
r(p) = [ dpe s{p) . (B.5) 



L is the common radius of AdS*, and the Sasaki-Einstein M5 in the solution dual to the 
superconformal theory, and is fixed by the number of D3 branes and the volume of the 
Sasaki-Einstein manifold. For T 1,1 : L 4 = ^-7rN c . 

The real integration constant C\ discriminates different classes of solutions. 

If Ci = then we recover the D3-branes solution (with nonzero c 2 , that can be fixed to 1) 
or its near- horizon AdS solution (with C2 = 0): 

e f = e 9 = e p = r (B.6) 
h = c 2 + ^. (B.7) 

If Ci > the solution describes iV c smeared D3-branes on the blown-up 4-cycle of the 
Calabi-Yau [3U[35|. Indeed, let us consider the change of radial coordinate: 

[l + ctf-*]- 1 = l- ^ = k{r) , (B.8) 



with r > b. If we further identify 
it follows that 



6 2 = cl /3 , (B.9) 

e 29 =r 2 (B.10) 

e 2 / = r 2 fc(r) (B.ll) 

e *V = , (B.12) 
so that the 6-dimensional metric, which is Calabi-Yau, is 

ds 2 6 = [k(r)]" l dr 2 + ^^(# + y cosOM 2 + % V (d6 2 + sin 2 6 t dtf) , (B.13) 

i=l,2 i=l,2 

that describes a deformation of the Calabi-Yau where a Kahler-Einstein 4-cycle is blown up 
at r = b. In order for the resolved Calabi-Yau to be smooth, an orbifolding along the U(l) 
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fiber parameterized by ip is usually needed. For the case of the deformation of the conifold 
the orbifold action is Z2, so that ip ranges from to 2n. The 10-dimensional metric of the 
solution is then: 

ds 2 10 = [h(r)]-V*dxL + Ur^dsl , (B.14) 



with the warp factoi 23 ' 



1, (r 2 -l) 3 1 /7T 2f 2 + l 

- log — 1 = arctan = — 

6 S f 6 -l v^V2 y/3 



r 

6 ■ (R15) 



The gauge theory dual to this local Kahler deformation of the Calabi-Yau cone is a deforma- 
tion of the superconformal theory due to the insertion of a VEV of a dimension 6 operator, 
which is a combination of the operators Tr(W a VV a ) 2 , W a being the gluino superfield [35J. 
The orbifold action is needed to have a dual field theory whose mesonic branch of the moduli 
space is (the symmetric product of N c copies of) the resolved Calabi-Yau. 

A similar analysis can be done for the solutions with c\ < 0, but in that case the 6- 
dimensional transverse space happens to have a curvature singularity and cannot be described 
as an algebraic variety. Therefore the supergravity solution is not expected to describe a 
dual supersymmetric gauge theory. 

Let us now study the unflavored limit of the general deformation of the KW solution of 
section [331 Recall that, in this case, the metric depends on three functions (/, g\ and #2) 
and the warp factor h. In terms of the variable p introduced in f l3.86f) the metric can be 
written as: 

ds 2 = h'^dxi 3 + ^ e 2f {dp 2 +l Y e 2 ^' 2f (de 2 + sm 2 e i d V 2 ) + \(d^+ V cosM^ 

6 ' 9 ^-^ 

^ i=l,2 i=l,2 

(B.16) 

The unflavored limit of the BPS system (13.851) amounts to taking C — 0. As in the previous 
case, the solution of section [3~5l is not valid in this limit (see eq. (I3.88P ) and one has to 
take C = in the system ( 13. 85ft and integrate it again following the same steps as in section 
13.51 The result can be written in terms of the function e 292 ^ p \ which is the solution of the 
following cubic equation: 

e 692 + ^a 2 e A92 = e 6p + c x , (B.17) 

where c\ is an integration constant. In terms of e 292 ^ the other functions of the ansatz can 
be written as: 

e 2 ^ = e 292 + a 2 , 



e 



2/ 



e 



6p 



e A 9 2 +a 2 e 2g 2 ' (B.18) 



h{p) = ~ Q I estate** + ° 2 



23 The additive integration constant in h is omitted in order to asymptote to AdS§ x X$ for large values 
of r. 
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where Q is given by (I3.84p and c<i is a new integration constant. Let us now perform the 
following change of radial variable 

2 

e 2 92 ( P ) = r _ _ (B.19) 
6 

Taking into account ( IB. 171) . the inverse relation between these two radial variables is: 

1 



6p 



216 



[ r e + 9a V - o 6 ] 



where we have redefined the constant C\ as: 

b 2 = 6(ci)3 . 

By using these relations, one can readily prove that: 

^91 . 



;b.2o) 



(B.21) 



- (r 2 + 6a 2 ) , 
6 



,2.f 



6 V ; 



where the function hz(r) is defined as follows: 

r 6 + 9a 2 r 4 - b 6 



K(r) = 



It is also easy to verify that: 



dp 



r 6 + 6a 2 r 4 
dr 



rn{r) 



;B.22) 



(B.23) 



(B.24) 



Using these results and redefining the warp factor as h{r) 2 — > h{r) 2/6, we get a metric of 
the form: 

(B.25) 



ds 2 = [h(r)] 2 dx\ z + [/i(r)] 2 ds 



with ds\ given by: 



ds 2 ; = [«(r)] X rfr 2 + ^- «(r) (# + ^ cos^d^) 2 + 



o=l,2 



+ - ( r 2 + 6a 2 ) ( d# 2 + sin 2 a dy? 2 ) + - r 2 ( d^ + sin 2 9 2 d<pl ) , (B.26) 
6 6 



while the warp factor h can be represented as: 



h(r) = -Q J 



rdr 



r 6 + 9a 2 r 4 - 6 6 



+ c 2 



(B.27) 



This is the solution with two Kahler deformations found in ref. [31]: the a constant param- 
eterizes global deformations, while the b parameter corresponds to local deformations. 
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C Appendix: Alternative Interpretation of the IR 
Regime 



Here we put an alternative description of the IR theory as deduced from supergravity, which 
honestly we could not discard. It mainly arises from the analysis of the Klebanov-Witten 
model at small values of the string coupling, and it is based on the non- validity of the orbifold 
relations fl2.63p - fl2.66p for all values of the parameters in the KW model, that was extensively 
pointed out in [13]. In the whole analysis that will follow, we will consider for clarity only 
the case of equal gauge couplings gi — g 2 = g. 

The curve of conformal points in the Klebanov-Witten model is obtained by requiring the 
anomalous dimension of the fields A, B to be 7a(<7jA) = —1/2, which assures (3 g = (3^ = 
(A is the dimensionless coupling from the quartic superpotential). The qualitative shape of 
the curve is depicted in Figure HI as well as some possible RG flows. The important feature 
is that there is a minimum value g* > that fixed points can have (due to the perturbative 
(5 g being negative, so that g = is an unstable IR point). One way to determine this 
curve of fixed points is to apply the a-maximization procedure originally spelled in [27] by 
using Lagrange multipliers enforcing the marginality constraints [IB], and then express the 
Lagrange multipliers in terms of the gauge and superpotential couplings. This computation 
for the Klebanov-Witten model was done in [49] o One can show that the curve of fixed 
points does not pass through the origin of the space of Lagrange multipliers, which is mapped 
into the origin of the space of couplings (free theory) . In a particular scheme the curve of fixed 
points is an arc of hyperbola with the major axis along A = 0. The exact shape of the curve is 
scheme-dependent, due to scheme-dependence of the relation between Lagrange multipliers 
and couplings: we choose a scheme in which the Lagrange multipliers are quadratic in the 
couplings. This choice fixes a conic section, and it is such a hyperbola because the one- 
loop anomalous dimensions of the chiral superfields get a negative contribution from gauge 
interactions and a positive contribution from superpotential interactions. The conclusion 
that the curve of conformal points does not pass through the origin of the space of coupling 
constants is physical. 

The family of KW SUGRA solutions describes the fixed curve. It is parameterized by 
that can take arbitrary values. For sufficiently large values of it, we can trust the orbifold 
formula: 2 

f- = e* for e+N c >l. (C.l) 

671 

The 't Hooft coupling g 2 N c is large (at least of order 1, so the theory is strongly coupled and 
the anomalous dimensions are of order 1) and the string frame curvature Rs ~ l/(e^N c ) is 
small. For smaller values e^N c < 1, ( jC II) cannot be correct: it would give small 't Hooft 
coupling while the gauge theory is always strongly coupled. The bottom end of the line 



4 We thank Sergio Bcnvcnuti for making us aware of this method and of the literature on the subject. 
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Figure 4: RG flow phase space for the Figure 5: Klebanov-Witten model with flavors. 

Klebanov-Witten model. The A-C flow has backreacting D7's in the A piece 

and then follows the KW line in the C piece; it cor- 
responds to Nf <C N c . The B flow is always far 
from the KW line, and corresponds to Nf > N c . 



corresponds to: 

{e*-0} <- {g = 9if ,\ = 0}, (C.2) 

and the SUGRA curvature is large even if the field theory is still strongly coupled. Anyway 
some quantities, for instance the quantum dimension of A, B, are protected and do not 
depend on the coupling, so they can be computed in SUGRA even for small values of e^N c . 

The supergravity solution of our system with D7-branes is in the IR quite similar to the 
KW geometry: the IR asymptotic background is AdS^ x T 1 ' 1 (with corrections), but with 
running dilaton. The field theory is thus deduced to be close to KW fixed line, but running 
along it as — > in the IR. Moreover, controls the gravitational backreaction of the 
D7-branes (as well as the gauge coupling), and as soon as e^Nf < 1 the branes behave as 
probes. In this regime, we expect the quantities computable from the background to be 
equal to the KW model ones: in particular 7^4 = —1/2. 

We can distinguish different regimes, starting from the UV to the IR. Depending on the 
values of N c and Nf they can be either well separated or not present at all. A section of 
the space of couplings and some RG flows are drawn in Figure El but one should include the 
third orthogonal direction h which is not plotted. 

• For 1 < we are in the Landau pole regime, and the dilaton (string coupling e^) is 
large. 

• For ^ < < 1 we are in a complicated piece of the flow, quite far from the KW fixed 
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line, as in the type A-B flows of Figure [51 In particular the D7-branes are backreacting. 
In this regime our SUGRA solution is perfectly behaved (as long as < e^). 

• For < < jg- (this regime exists for Nf < N c ) we are in a region with almost 
probe D7-braneo, so we are close to the KW line, but with large 't Hooft coupling, 
so we can trust (IC.ll) . We can expect the energy /radius relation to be quite similar to 
the conformal one, thus we can compute the gauge /3-function and deduce the flavor 
anomalous dimensions 7q. Apart from corrections, we get: 

11 13 

Ja — ~2 Ra ~2 1q ~ 4 Rq ~ I ' (°' 3 ) 

The R-symmetry is classically preserved but anomalous as in supergravity. The various 
/5-functions are computed to be 

We want to stress that this regime in not conformal, and in fact the theory flows along 
the KW fixed line, as in the type C flow of Figure [51 The smaller is Nf/N c , the longer 
is this piece of the flow. For Nf > N c this regime does not exist, and the theory follows 
the type B flows of Figure [51 

• For < Min(^, -^) we are close to the end of the KW fixed line, and the gauge 
coupling is close to g*. Again the D7's are almost probes. The string frame curvature 
is large, as in the KW model at small g s N c . Since the gauge coupling cannot go below 
g*, its /5-function vanishes even if the string coupling continues flowing to zero. We get 
in field theory: 

IQ - 1 Rq * \ (C5) 
f$n = \h. (C.6) 

All the flows accumulate at the point {g — g*, X — 0} of Figure El but the theory is not 
conformal. In fact the coupling h always flows to smaller values, and the theory moves 
"orthogonal" to the figure. For this reason 7q and Rq do not satisfy the relation of 
superconformal theories. 

• The end of the flow is the superconformal point with h = (and g = g#), which should 
correspond to = and cannot be described by supergravity. Without the cubic 
superpotential one can construct a new anomaly free R-symmetry with Rq = 1, by 

25 The dual in field theory of the D7's being probes is that graphs with flavors in the loops are suppressed 
with respect to gauge fields in the loops, since Nf < N c . 



1 p 1 

P 9 ^0 ft 
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combining the previous one (Rq = 3/4) with the anomalous axial symmetry which 
gives charge 1/4 to every flavor. This satisfies known theorems on superconformal 
theories. Moreover, the fact that h — > in the far infrared realizes in field theory the 
incapability of resolving the D7 separation at small energies, and the flavor symmetry 
S(U(N f ) x U(Nf)) is restored. 

Note that when Nf > N c and the D7-branes are probes (this is the regime < < jj- 
and g = g*) one could think hard to see in field theory a suppression of graphs with flavors 
in the loops, with respect to gauge fields in the loops. Consider the gauge propagator at 
1-loop with flavors (Figure EJ). It is of order g*Nf, not suppressed with respect to the graph 
with gauge fields in the loop of order glN c . But if we sum all the loops with flavors, we must 
obtain the flavor contribution to the /5-function, which for g ~ and so 7q ~ 1 is indeed 
very small. 




7 = 1 7Q = i FLOW Pole 



N c N f 



Figure 6: Flavor 1-loop correction to the Figure 7: Regimes of KW with flavors for 
gauge propagator. Nf < N c . 

A summary of the phase space for Nf < N c is in Figure [71 The computation in [2T] is 
valid in the region < < 4- of the phase space. 
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